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COMMUTATIVE POLYNOMIALS 


By H. D. BLOCK and H. P. THIELMAN (Ames, Iowa) 


[Received 13 March 1950; in revised form 18 December 1950] 


We define the functional product of two functions f(x) and g(x) as f[g(x) | 
and denote it by fg. The ordinary product of two functions f(x) and g(x) 
will be designated by f(x).g(2). We use g” to denote the nth functional 
iterate of g, while [g(x)|" denotes the ordinary nth power of g(x). The 
function f(x) = x will be called the identity.* If fg = gf = x, we say 
that g is the inverse of f, and write g = f-!. We write g,, with the 
subscript denoting that g, is a polynomial of proper degree k. 

Two functions f and g are said to commute with each other if fg = gf. 
An entire set of commutative polynomials is a set of polynomials which 
contains at least one of each positive degree, and is such that any two 
polynomials of the set commute with each other. We confine ourselves 
to polynomials with real coefficients. It is the object of this paper to 
give all entire sets of commutative polynomials and some of their proper- 
ties. Some of the results given here can be derived from certain results 
of Ritt,+ but they were not explicitly stated in his papers. The methods 
used in those earlier papers involve, however, more advanced work in 
analysis than ours. To establish these results by entirely elementary 
means we give some preliminary definitions and a lemma. 

If = ax*?+ then we define A(h,) = B?—4ay. 

Lemma. If h,(x) commutes with hg(x) then A(h,) = 28 or 28-+-8. 

This lemma is easily established by equating coefficients of like powers 
of x in h,h,(x) = 

If A(z) = Ax+ B(A + 0), then A-'fA is called the linear transform of 
fby A. It is readily seen that, if f and g commute, then so do A-"fA, A-1gA. 
From the fact that the set p,(2) = x” (n = 1, 2, 3,...) is an entire com- 
mutative set, and /,,(”) = cosn(arccos 2) is also an entire commutative 
set we can define the entire sets of commutative polynomials: 

(1) = = {(Az+B)"— BY A, 

(II) = cosnarccos A(x) = A-!cos n{arccos(Ax+ B)— B} 

(A #0; = I, 8, 3,...). 


* We shall use the notation and terminology that A. G. Walker uses in the 
paper ‘Commutative functions’, Quart. J. of Math. (Oxford) 17 (1946), 65-92. 
{ Ritt, J. F., Trans. American Math. Soc. 23 (1922), 51-66 ; 25 (1923), 399-448. 


Quart. J. Math. Oxford (2), 2 (1951), 241-3. 
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242 H. D. BLOCK AND H. P. THIELMAN 
THEOREM. The only entire sets of commutative polynomials are the sets 
(I) and (II). 


Proof. It is readily verified that A(h,) = 28 if and only if we can find 
a A such that h, = Similarly A(h,) = 28+ 8 if and only if we can 
find a A such that h, = A-1cos 2(arccosA). In view of the lemma, if 
h,,(k > 3) belongs to an entire set of commutative polynomials, then 
h, commutes either with P, or with T,. We shall now show that the 
first possibility implies while the second implies h; = thus 
completing the proof of the theorem. 

If h, commutes with P,, then H;, = Ah, A~! commutes with x?. Thus 


= (1) 


for all positive integers n. Equation (1) shows that every zero of H,,(x?") 
is also a zero of H,(x). Thus, if x, = re“ (r 4 0,0 ¢ 0) is a zero of H,,(x), 
so is r/2ngi/2n This gives an infinite number of distinct zeros for H,,(z), 
which is impossible. Therefore every zero of H,,(x) is the number 0, and 
H,,(«) = ax*, Equation (1) then shows that a = 1. Finally 


h, = = P,.. 
If h, commutes with T,, then H;, = Ah, A~-! commutes with 
cos 2(arccos 2) = 2a?—1; 


i.e. H,(2a2—1) = 2[H,(«)?—1. (2) 


Differentiating both sides of (2) gives 

a.Hi(2a2—1) = H,(x).Hi(a). (3) 
Now let F(x) = P+ (@) (4) 
Making use of (2), (3), (4) we find 

F(2x?—1) = 4[H,(x)]?. F(a). (5) 


Equation (4) shows that F(x) is of degree less than 2k. But, if 
F(x) 4 0, then equation (5) shows that the degree of F(x) is 2k. Thus 
F(x) = 0, and (4) becomes a differential equation in H,(x), which com- 
bined with (2) shows that H,(1) = 1. The only polynomial solution of 
the differential equation satisfying this condition and (2) is 


H,,(x) = cos k(arccos x), 


and thus = T,,. 
It might be of interest to note that the polynomials of the set (II) all 
have the same fixed point 6 = (1—B)/A; ie. 7,(b) = 6. If we let 


COMMUTATIVE POLYNOMIALS 243 


a = 2A, then these polynomials, expanded about the fixed point 5, 
take the form 


T,,(x) = 20-1 cos n{arccos }(ax—ab+ 2)}+b—2a-1 


This last equality can be verified by mathematical induction assuming 
its truth for n—1 and n, and showing that it then holds for n+1. Since 
TT, = TT, = Tn; We obtain the identity 


k=1 j= 


akA(nm+k—1)! 
(nm—k)!(2k—1)!k 


which might be quite difficult to establish some other way. 

We might just add that it follows from the results of Ritt* combined 
with ours that, if two non-linear polynomials commute with each other, 
then (except for certain ‘trivial cases’) they both belong either to set (I) 
or to set (II). The ‘trivial cases’ to be excluded consist of the following 
obviously commutative sets: (i) {f"(x)} (n = 1, 2,...); (ii) if g commutes 
with / and if g is an odd function, then g commutes with —)h; (iii) linear 
transforms of the cases given in (i) and (ii). 


* Loc. cit. 


SERIES EXPANSIONS OF CONTINUOUS GROUPS 
By D. D. KOSAMBI (Bombay) 


[Received 17 Ap.il 1950; in revised form 2 December 1950] 


In this note, I consider continuous groups of transformations possessing 
a unique formal expansion in powers of the parameter, but regardless 
of the question of convergence and essentially in the neighbourhood of 
the identity. 


1. Suppose a vector space V of elements ¢ to be defined over a con- 
stant field C, which I shall assume for simplicity to be that of all real 
numbers. By definition, a continuous one-parameter group is the map- 
ping ¢— 7(¢,t) which maps V into itself for each value of the real 
parameter ¢ in some open ft-interval containing the point t= 0. It 
follows that this t-interval is independent of the particular element ¢, 
though the transform 7’ may be defined for other values of t outside the 
interval for some sub-spaces of V. Though the treatment that follows 
is essentially local, this assumption is a step towards consideration of 
the whole group, being of more than local character. However, in the 
formal power-series expansion, which is the main tool of analysis (as in 
the classical work of Liet and his followers) the interval of convergence 
will not be discussed, thus restricting the treatment again to local 
phenomena, namely the ‘infinitesimal’ part of the group. Finally, we 
consider only the case where, for suitably restricted ¢, no two t-values 
give the same mapping. This excludes the periodic case as well as the 
trivial case where the entire group is the identity. The usual group 
properties are now demanded, as follows. 


(i) The group law. This reads 


T(T(¢,t), 8) = T(¢,$(s,t)). 
Or, abbreviating 7'(¢,t) as 7'(t), 
T(s)T(t) = T((s,t)). 

That is, two successive transformations represented by values t, s of the 
parameter being performed in that particular order, they have a combi- 
nation which is another transformation of the same sort with parameter 
Y(s,t). Hypothesis: (s,t) is a continuous function of the two variables 

t Summarized in 8S. Lie and F. Engel, Theorie der Transformationsgruppen, 
vols. 1-2 (Leipzig 1888) ; vol. 3, 1893. 

Quart. J. Math. Oxford (2), 2 (1951), 244-57. 
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s, t, uniformly approximable in some neighbourhood of (0,0) by power 
series themselves satisfying the special properties which may be de- 
manded of (s,t) in the sequel. We proceed to investigate the nature 
of the function , again regardless of what happens when ¢ is changed; 
from the present point of view, the relevance of this latter question 
arises only when the change is due to another group, which may be 
represented by the introduction of another parameter independent of t. 
The invariance of sub-spaces of V under the group or some transforma- 
tions thereof is also a question of which no cognizance is taken in the 
treatment. 


(ii) The identity I. This we take as fixed once for all as at t = 0, so 
that 7'(¢,0) = ¢ for allde V. By known results in elementary group 
theory, the identity is both a right and a left identity, whence 


(0, t) = P(t, 0) = t. 
The approximating expansions to % must therefore have the form 
b(s,t) ~ s+t+ast/2!+ (bs*t+-cst?)/3!+.... (1.1) 


(iii) The inverse. To every t there corresponds one and only one real 
value of the parameter s such that ¢(s,t) = 0. 

The real solution of y(s,t) = 0 for s as a function of t is therefore to 
be single-valued. But, again, it is known that the right and the left 
inverse are identical, whence this solution coincides with that of 
u(t,s) = 0. If the solution in question be s = f(t), it follows that 
t = f(s) and f(f(x)) = & for all real x. By previous considerations, the 
identity being its own inverse, f(z) must be continuous at x = 0, with 
f(0) = 0. Assuming formal expansibility of f(x) in power series shows 
that the only solutions are then f(x) = +a/(cx1); which can be re- 
duced by a schlicht analytic transformation to f(z) = +2, and we may 
further assume this to have been done. The plus sign implies that the 
transformation is its own inverse, whence 7? = J for every permissible 
value of the parameter, so that (t,t) = 0. For any two transformations 
S, T of the group, 

(7S)(TS) = TSTS = I; 
but 7 = T-1, S = S-1, whence 


T-*§ = = 8T-, 


the group being Abelian. Therefore = s). The only approxi- 
mating function satisfying all these conditions on ¢ is zero; whence the 
whole group is the identity, a trivial case that has already been excluded. 


‘ 
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In its standard form, with the identity at t = 0, T(¢,t) has the inverse 
T(¢, —t). 
The expansions for ¥(s,t) now reduce to 
t) ~ (s+t){1+b, st+(b, s*t+-c, st?)/2!+...}. (1.2) 
(iv) The associative law. This has already been utilized in the pre- 
ceding section. Our statement is 
T(s{T(t)T (u)} = {7 (s)T()}T(u) = T(s)T(t)T (uw). 
This has as consequences: 
The second is obtained by taking 7(u) = T-1(t)T-(s) in the first. 
Further, by taking t+-u = 0, we obtain the two relations 
t), —t} = o{—t, 8)} = s. (1.4) 


Noting that, for any fixed ¢, #(s,t) and Y(t, s) run through a continuous 
range of parameter values which we may call w, w’ respectively, we 


H(w,t) = Yt, = —W(—t, —w). 
These being again identities, one may compare the two sets of approxi- 
mating expansions to reach the conclusion that w = w’, i.e. 
t) = P(t, 8). 
Thus we have shown that 


every one-parameter continuous group is Abelian, and the combination 
function x(s,t) is of the form (s+t)~’, 
where —t) = w'(t,s) = 

2. Rather than discuss the functional equations for »(s,t) I proceed 
to further standardization by using the properties of the vector space, 


and the main hypothesis: The general transformation T(¢,t) is uniquely 
expansible as a formal power series in t 


T($,t) = $+tK, (2.1) 


This series representation is certainly valid at t = 0, giving the iden- 
tity, as required. The case of finite polynomial expansions presents no 
difficulty. Where the vector space has a norm, or some other means of 
defining convergence, the expansion has the usual meaning. But even 
when the series does not converge, or in the case of vector spaces for 
which no norm is defined so that convergence cannot be discussed at all, 
the formal power-series representation may still be used correctly, by 
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analogy with modern procedure in algebra and analysis,t under the 
following restrictions: 


(a) Every expansion for a 7(¢,t) is to be unique, even when the 
mapping depends upon a finite number (greater than unity) of para- 
meters, and the expansion is in multiple power series. (5) If, starting 
with a given group and its expansion, we substitute for the parameter 
an expansible function of one or more variables, the resulting formal 
expansion in powers is still to be the unique representation. (c) Every 
operation that is to be performed upon an element of the group may 
also be performed (as for a sum) upon its series representation, in the 
sense that the resulting series then represents the group element conse- 
quent upon the said operation. 


In essence, we retain all the properties of the series as a device for 
manipulation when associating a set of operators K,, Kg,... with a given 
transformation 7'(t). However, the ability to say whether or not an 
arbitrarily given series expansion represents some element of the group, 
or even of the vector space, has been lost. 

The coefficients of (2.1) depend upon ¢ alone, not t, by hypothesis; 
the restriction that the representation have the precise meaning of 
equality in the case of finite expansions or convergence makes it neces- 
sary that the K,, ¢ be also elements of V for all n, seeing that 7'(¢, t) is 
to be a mapping of V into itself. We may exclude any representations 
where the coefficients of the power series are not elements of ¢, though 
the possibility is of interest. Therefore, the K; are operators mapping V 
into itself. From the group law we get the following identity: 


T(s)T(t) = $+tK, 
+-8K,(6-+tK, +... 
= $-+(s,)K, (2.2) 


It follows that the group product 7'(s)7(t) is expressible as a double 
power series, and the coefficients on both sides are equal (by our assump- 
tions above) term by term for each product s?t?. We may assume (f(s, t) 
itself to have a formal power-series expansion without loss of generality, 
and in fact have virtually done so in the previous section. 


+ As, for example, in asymptotic series, which need not converge anywhere but 
are useful both for numerical calculation and for representation of functions 
allowing formally algebraic operations, differentiation, etc., cf. K. “nopp, Theory 
and Application of Infinite Series (London, 1928) §§ 300-2. 


j 
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The operators K themselves define a vector space over C with reference 
to V. That is K;¢, aK;¢, (aK;+-0K;)¢ are all in V if V, anda, be C. 
Thus, any finite linear combination of the operators is defined, as well 
as the product K;K;, by means of K;,(K;¢); the algebra of operators is 
power associative. For brevity, the tensor summation convention is 
introduced here so that an index repeated in subscript and superscript 
is to be summed over all possible values; ¢”K, will be the whole formal 
infinite series from values 1 to 00 of the index. The equality (2.2) says 
¢/m!) = ¢, 
where the a?; must be independent of ¢t and ¢, but in C. For such fixed 
constant coefficients from C, the powers of t over any t-interval whatso- 
ever about the origin are linearly independent. For special non-zero 
values, say t,, ty,..., the set of operators ?’K,/r! (which exist by (2.1)) 
would have the same linear dimension as the original set K,, Kg,.... 
Choosing a suitable basis from these, or even directly, we see that each 
operator K is linear when applied to the operator-space of K,, Kg,.... 
The crux of the argument derives from assumption (c) above. When 
the identity operator J is included, we have at least a restricted linearity 

K,(1+aK;) = K;+aK,K; (ae C). 

Thus we have shown that 


The operators K in the expansion of a one-parameter group in formal 
power series have the property of restricted linearity over the operator-space 
defined with reference to the original vector space in the sense that 


and full linearity over the operator-space excluding the identity. 

As a corollary, it follows that it suffices for the K; to be all linear in 
the ordinary sense; also that complete linearity follows if the union of 
all sub-spaces into which V is mapped by the K; exclusive of the identity 
coincides with V. 

This still does not exhaust the content of (2.2), which also determines 
the actual algebra of operators, by giving explicitly each product K; K;. 
First rewrite (2.2) as 


b+(s+t)K, 6+{(s+t)?K, 6+ 
K,—K;)}/3!+... 


= (2.3) 
This leads to 


Ki-K,=0, XK,K,=K,K,=K,  6,K,+3K,—3K? = 0, etc., 
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where 0 is the null-operator and the coefficients 5,,... are precisely those 
which occur in the expansion of ¢(s,t) in (1.2). 


K,, is expressible as a polynomial of degree n in K, with the coefficient 
of the highest power equal to unity and no ‘constant’ term. 


Inasmuch as the definition with reference to V of powers of K gave 
power-associativity, we could have proved our previous result again, 
namely that every expansible one-parameter group is Abelian, whence 


= w(t, 8). 


3. A further canonical reduction is now made by replacing ¢ with a 
suitably expansible function g(t). This means replacing K,, (n > 1) by 
linear combinations of the K, (r < n) still keeping K, unchanged, intro- 
ducing no constant term, and retaining unity as the coefficient of K,,, 
provided that 


g(t) = t+cgt?/2!+-c, #3/3!+.... (3.1) 
This will give, according to our previous hypotheses, 
T(¢,9(t)) = 6+tK, Ky)p/2!+.... (3.2) 
The plan is now to show that, by suitable choice of the function g(t), 
the operators become K,, K?, K3,.... For this we must have 


Kj? = Ki+c, K, and hence c, = 0, excluding the trivial K, = 0; 

K? = K3—b, K,+c, K, and hence c, = and so on. 
It can be verified that the coefficients c are uniquely determined as 
linear combinations of such coefficients as survive in the expansion of 


%(s,t). But after the transformation, we have the new % = s+, as can 
be seen directly. 


There exists a formally expansible transformation g(t) of the parameter t 
such that the canonical form becomes 


T(d,t) = = (exptL)d, 


L being an operator with restricted linearity over the identity and positive 
integral powers of L, defined with reference to V and C. 


As a corollary, we get another classical result: every one-parameter 
transformation group becomes a translation in the (standard) parameter. 
This is a consequence of ¢'(s,t) = s+t after standardization. The exis- 
tence of a formal series for g(t) such that 


%(9(s), g(t)) = g(s+t) 


could have been shown directly at the end of the previous section. 


: 
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I consider two examples: 


(i) Let the group be represented by the transformations 
w = (z—t)/(1—tz) #1). (3.3) 


This is familiar as that linear fractional transformation of the complex 
plane which carries the point z = ¢ (under our restrictions on the real 
axis) into the centre of the unit circle, while carrying that circle into 
itself. These transformations constitute a continuous group with the 
identity at the origin t = 0, the composition function being 


t) = (s+t)/(1+st), 
which can be reduced to the form s’+-t’ by means of the familiar trans- 
formation t = tanht’. The formal series expansion of (3.3) is 


w = (3.4) 


At first sight the operator Z would seem to be that which replaces z by 
z*—1, but it is then seen that even the third term is not correctly 
expressed in the desired form. The fact is that we have not yet con- 
sidered the vector space, which is that of analytic functions of z, and 
not just the real axis in the z-plane. The correct operator for the 
canonical expansion is (z?—1)d/dz, and the direct expansion is imme- 
diately to be expressed in terms of polynomials of this. It will be shown 
later on that the general L corresponds to a differentiation. 


(ii) The vector space consists of integrable functions ¢(x) of a real 
variable x on some fixed set of points and for some fixed process of 
integration. By ‘integration’ is meant a mapping of the class of 
integrable functions into the set of real numbers, the operation being 
completely additive over disjoint subsets of the fundamental set; being 
linear distributive over C; and non-negative if the integrand is non- 
negative, over any such subset. If the characteristic functions of a 
class of subsets be included in the class of integrable functions, the 
subsets will be called measurable, and the evaluation of the integral may 
then proceed as usual (with the further requirement of separability, etc.) 
using this measure. In addition, we can follow the procedure given by 
L. Schwartz in generalizing the concept of ‘function’ by assigning to 
an abstract function f a class of mappings i f¢ defined for the whole set 
of integrable functions ¢. 

In this abstract and general sense, the operator L is to be defined by 
means of a kernel Lé = | K(x, y)¢(y) dy, being therefore linear, and a 
mapping which depends upon 2, i.e. a mapping of the function-space 
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into itself. Defining L*, L%,... by means of the iterated kernels 
LK = | K(x,w)K(w,y)dw, = LK-1, 
we have the group given by 
T($,t) = | K(x, yy) dy +l K*(x, y)$ly) dy/2!+.... 


Ifthe kernel K has a discrete spectrum and the corresponding expansion 
in (orthonormal) eigenfunctions K = > A; f,(x) f(y), we have 
T(¢,t) = $*(x)+ a,(expt; f(x), 

where the a; are Fourier coefficients of ¢ with respect to the functions f,, 
¢* being that part of ¢ which lies in the section of Hilbert space ortho- 
gonal to the portion swept out by the set f;,, being equivalent to zero 
(with respect to the process of integration) if the set f; is complete. 

We may also use these ideas to obtain a different type of group expan- 
sion, arbitrarily defining 

T'(¢, t) o+tLd, T(s)T(t) o+(s+t) Ld. 

This means defining the operator algebra without reference to V, so that 
= 0 for all whence L(L¢d) ~ in general. In such definitions 
the advantages of power associativity will have to be obtained by special 
devices. Here, the group is in effect a generalized translation in the 
vector space; all previous developments can be restored if there should 
exist a class of elements » in V such that Ld € {u}, Lu = 0; the transla- 
tion group is then simply written as 7'(¢,t) = ¢+tp, uw being a particular 
element of the class. 


4. It is now easy to develop a series of results which give the actual 
determination of the canonical operator. The first is an immediate con- 
sequence of the steps used in reduction to the standard form: 


The operator L which gives the formal canonical expansion — 
T(¢,t) = (exptL)¢ 
is to be obtained by reducing T(¢,t)—¢ modulot whenever the identity is at 
t= 0. In fact, we saw that after taking the identity at ¢ = 0 the further 
reductions could not affect K, = L. This operator L will be called the 
infinitesimal transformation of the group, seeing that the canonical form 
is its exponential mapping, at least in the formal sense. 

Suppose that ” essential parameters are involved instead of the single 
parameter to which the discussion has hitherto been restricted. We 
standardize as before by taking the group-identity as at the origin in 
the parameter n-space. Repeating the previous arguments, and taking 


4 
| 
as 


252 D. D. KOSAMBI 


account of the fact that the group is also a one-parameter group in each 
of the arguments, we find at once that there must exist » operators 
L,, L,,..., L,,. These must be linearly independent over C, or the total 
number of parameters will be reducible below the n which we have taken 
as essential, thus contradicting the hypothesis. It is clear that a canoni- 
cal form exp(t, L,+t, L,+...+t,, L,,) does not exist except in the Abelian 
case. But, if a transformation of the group is built up by performing the 
transformations associated with ¢,, t,,... in some fixed order, and then 
we replace all the parameters by expansible functions of a single para- 
meter ¢, we should again have a one-parameter group which can be 
standardized. It is seen at once that the infinitesimal transformation 
of this simplified group will necessarily have an operator 
a, L,+a,L,+...+a, L,, 

where ¢t; = a;t+... are the substitutions. On comparing the original 
multiple expansion with the new one, it will be seen that equality is 
possible if and only if L;L,—L,L, is a linear combination of the 
operators with coefficients in C, for every pair of indices i + j. 

A finite set of operators (defined with respect to V, with the property of 
restricted linearity) linearly independent over C and forming a Lie algebra 
is defined by every expansible multiparametric group. 

The converse cannot even be approached by our present methods, for 
the formal expansions may be set up and shown to satisfy all the group 
conditions provided the operators do form such a Lie algebra; but there 
is no method by which the expansions can be shown to represent any- 
thing in V. To show that a Lie group exists for every such Lie algebra 
will obviously need a different approach. 

The operator Z has some of the properties of differentiation besides 
restricted linearity: 

Suppose a product ¢.x defined as a mapping of the Cartesian product 
V x V into some vector space P over C, with the properties 


(ag, = € P, 
(4x, +bx2) = a,beC, (¢,x)EVXV. 
Then the group in V with infinitesimal transformation L induces another 
in P with infinitesimal transformation 
L(g. x) = (L4).x+4.(Lx), 
provided that the latter group is also formally expansible. 


For proof, it has merely to be noted that the group in P is the mapping 
¢.x > T(¢,t).T (x,t). The unique expansibility of this, by hypothesis, 


a 
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and the definition of the product allow us to multiply the corresponding 


power SNS = 
which has only to be reduced modulo ¢ for the new infinitesimal trans- 
formation. Actually, we could extend this to two separate standardized 
groups in a more general product space 
x) = x), 

with obvious modifications in the definition. We do not necessarily have 
la = 0, which is so essential a property of ordinary differentiation. 

Let V be the real line with coordinate x, and f(x) an analytic mapping 
of V into itself. Then we know that 

fle+t) = flw)+tf' (a) 

The derivative is the infinitesimal transformation which generates the group 
induced upon the vector space of analytic mappings when the base-space 
undergoes a translation. 


This is immediately extensible to n dimensions, and can be used also 
as our definition of the derivative. From this and the definition of an 
infinitesimal transformation, the further result is obtained: 


Given a vector space of real analytic mappings f(x',...,x") of R,, into R,; 
the operator generating the group induced upon f by a one-parameter trans- 
formation group x‘ > T(x',t) with the infinitesimal transformation 

Lat = u*(z',..., 2") 

is Lf = = u'f,. 

(The very last portion merely explains a notation standard in differential 
geometry, the tensor summation convention being again used while the 
comma in subscript denotes ordinary partial differentiation.) The result 
follows from the fact that we may regard the transformation upon each 
as a translation = Use of the previous theorems 
plus the reduction modulo ¢ applied to T(f,t)—f yield the formula. This 
and the properties of the product operation will be useful in deriving 
what is essentially the Lie operator in abstract differential geometry. 


5. The foregoing may now be used to evaluate L for special types of 
vector spaces. Let us first suppose V to have been mapped onto an 
n-dimensional analytic manifold in the sense of Whitney—Chevalley.t 


+ C. Chevalley, Theory of Lie Groups, vol. 2 (Princeton 1946) Chapter III; 
H. Whitney, Annals of Math. (2) 37 (1936), 645-80. This would restrict our vector 
space by means of topological considerations, which naturally lead to the study 
of the whole group. 


‘ 
L 
by 
| 
| 
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The last result above says that, for the analytic functions on that mani- 
fold, Lf = u’f, is an analytic function obtained by a linear operation. 
This means that w‘(x) is a contravariant vector over the manifold, i.e. an 
element of the space of all linear mappings of the partial derivatives f; 
into the domain of real numbers (more strictly, into the analytic mani- 
fold). Then f; is itself an element of the dual space, a covariant vector. 
Tensors of any rank, covariant as well as contravariant, can now be 
defined as the Kronecker products of a suitable number of spaces of each 
kind. The problem is to deduce the explicit formula for L in the general 
vector space which tensors of each fixed rank constitute. From the very 
definition, the tensors may be treated as symbolic products of a number 
of ideal vector elements equal in number to the rank. Using the theorem 
for the L-operation in a product space we have only to derive the actual 
operators for covariant and contravariant vectors. From the fore- 


going, = L(Lf) = wu'f,) 


But the term within the brackets is itself a scalar product to which the 
product theorem applies; so that 


L*f (Lu’) f,+u'L\ f, fp 


Now it may be verified as an elementary exercise in tensor analysis that 
neither of the two terms in the final expression is a scalar. Moreover, 
it can be rewritten as u’( f,, v4.+-u?f,,). The bracketed expression being 
a covariant vector can be identified only with Lf, We have thus 
proved that 


The L-operator has (besides linearity) the properties Lu’ = 0 where 
u'(x},...,2") is the contravariant vector that defines L; and 


Tr 9 
Ln, = 
for any covariant vector v;. 


The second part may be seen directly from the calculation for L?f or 
deduced from application of LZ to the scalar product u’v,, wu being the 
defining vector so that Lu‘ = 0. For the general contravariant vector 
w', operate once again on w'v,, where we use the convention of an upper 
index for contravariance, and the lower for covariance. This gives us 
Iv = w'w',—w'u',, which is precisely the Poisson bracket fu, w}. Now 
that L is known for both sorts of vectors, we use the product theorem 
to get: 
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The general L-operator for a tensor Ti: is given byt 

= Ti +...—uil, Ti. (5.1) 
The last term is for tensors of weight p (i.e. tensor densities), which are 
multiplied by the pth power of the Jacobian of the transformation, on 
change of point coordinates over a portion of the manifold. The classical 
tensor analysis deals solely with transformation laws for such coordinate 
transformations in defining tensors. From that point of view, the general 
formula for LT~ is merely the difference of the changes due to the infinitesi- 
mal transformation regarded first as a change of variable and then as a 
change of coordinates; this may be taken as the basic definition of the L- 

operator as applied to any geometric object, tensorial or not. 
By a geometric object is meant, in classical tensor analysis, any func- 
tion of the n variables which has a law of transformation also given when 
a change of coordinates is effected. For tensors, this is linear in each 


component: Ti = T* 


Here, the transformation #* = f*(z!,...,2”) is non-singular, with an in- 
verse, and the necessary differentiability is always assumed. The trans- 
formation coefficients satisfy = 8, and the contraction 
3; = n, the dimension of the space. There also exist geometric objects 
which have non-linear transformation laws, and may involve higher 
derivatives than the first. Such are, for example, the Christoffel symbols 
of Riemann geometry, included in their immediate generalizations, the 
coefficients of symmetric affine connexion I'},. These are used to define 
covariant differentiation; this, from our point of view, may be described 
briefly (if it exists) as an operator of partial differentiation which defines 
the tangent space for each tensor manifold (defined upon the original 
analytic manifold) with the same set of differentials. Substitution of 
covariant for ordinary partial differentiation in (5.1) leaves the formula 


+ The first known case of this operator being used in differential geometry is 
presumably in the famous ‘equations of Killing’ equivalent to Lg;; = 0, Crelle’s 
Jour. 109 (1892) 167. This, incidentally, assumes, as usual, that L#' or rather 
L dx‘ = 0, though the validity of the full series expansion is nowhere considered. 
The first explicit statement of (5.1) seems to be in J. A. Schouten and E. R. van 
Kampen, ‘Beitrage zur Theorie der Deformation’, Wydawnictwo Redakcji Prac 
Mat.-Fiz. (Warszawa 1933) 4-5, where LT}, is also given, but without recognizing 
its connexion with collineation groups or the series expansion. The L-operator 
for tensors of rank zero but weight not zero may be seen in the preservation of 
volume for the motion of an incompressible fluid, where the lines of flow are 
really the trajectories of a one-parameter group ; for the infinite series expansions, 
ef. D. D. Kosambi, Math. Student (1934), 94-100 and M. Shabbar, J. Indian Math. 
Soc. 6 (1943), 186-91. 


Fj 
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for LT; unchanged. The advantage of extending the L-operator to such 
general geometric objects is that it saves a great deal of labour in the 
study of groups associated with various types of differential geometry. 
In particular, LT, = 0 gives the (infinitesimal) collineation group for a 
Riemannian or an affinely connected path-space.+ These equations are 
short-hand for a system of second-order partial differential equations 
whose vector solutions are precisely the u‘ that define the group- 
operators. Moreover, the compatibility conditions for these partial 
differential equations are invariably reducible to the preservation under L 
of the differential invariants of the space. For LY‘, = 0 we get the con- 
ditions of integrability as the first-order system LRi,, = 0, where R is 
the curvature tensor; for the projective collineations, LWi,, = 0, W 
being the Weyl tensor, and so on, as may be seen either by reducing 
the formulae given by various workers in this field, or from first con- 
siderations as above. 

Such groups of collineations are in fact the solutions of the equations 
of variation.t For example, the second-order system 


= 0, 
where ¢ is an absolute parameter and the dots denote differentiation, 


includes both the geodesics of the Riemann space and the affine paths. 
Their equations of variation are 


dat = 0, 
where the dots on wu‘ indicate ‘total’ differentiation along the paths. 
If we look for the solutions u‘(x, #,t) we get a Lie ring for the extended 


+ L. P. Eisenhart, Riemannian Geometry (Princeton 1926), § 68; it may be 
pointed out that the linear equations in &,, just following § 68.14 are identically 
satisfied in virtue of known properties of the curvature tensor. For affine and 
projective collineations, Eisenhart, Non-Riemannian Geometry, §§ 46-9; equa- 
tions (47.6) are reducible to LWi,, = 0. 

t D. D. Kosambi, Proc. Nat. Acad. Sci., U.S.A., 35 (1949), 389-94. This shows 
that the full mechanism of tensor analysis usually employed to prove the group 
property of the solutions is not at all necessary. For ‘collineations’ in path-spaces 
defined by second-order partial differential equations (as a space of ‘K-spreads’) 
E. T. Davies, J. Lond. Math. Soc. 18 (1943), 100-7; Buchin Su, Trans. American 
Math. Soc. 61 (1947), 495-507 ; Science Record (Academia Sinica), 2 (1947), 11-19. 
For the tensor analysis of ordinary differential equations of any order and partial 
differential equations of the second order, D. D. Kosambi, Quart. J. of Math. 
(Oxford), 6 (1935), 1-12; 7 (1936), 97-104; 19 (1948), 204-19; V. Seetharaman, 
Proc. London Math. Soc. 2 (1938), 64-87; Proc. Indian Acad. Sci. 21 (1945), 
211-17. These develop the differential geometry without K-spreads or allied 
devices, and list the basis for differential operators as well as differential invariants. 
The problem of passing from local to the over-all topological invariants even for 
such specialized fibre-spaces remains unsolved. 
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operators L = wu’ 0/é%"+-u" é/éx". Solutions u‘(x) give the Lie algebra, 
and in the analytic case the Lie group preserving the base space. It may 
be noted that the general space is a type of fibre-space, with the ‘velocity’ 
vectors * constituting the fibre-bundle attached to each point of the 
base-space La’ = 0; along the paths we have a ‘total’ differentiation 
djdt = —of +0/ét, 
(which has to be used in the geometry and in L above); the equations of 
variation give precisely those extended operators L which commute with 
d/dt. These ideas and procedures are immediately extensible to higher- 
order ordinary differential equations, as well as to systems of partial 
differential equations. Our methods give the Lie algebra associated with 
a group, but the whole group only in the analytic case—though it seems 
obvious that the singularities of the differential equations must be re- 
lated to those of the group expansions. Almost any fairly complicated 
‘space’ admits only the identity, so that classification of systems of 
differential equations by means of the groups admitted is not definitive. 
Nevertheless, our methods may now be used in abstract differential 
geometry, even to carrying over the notation. The fundamental prob- 
lem of defining tensors in abstract spaces, differentials, and path- 
equations of various sorts has been dealt with by Michal} and his 
collaborators. The groups preserving such spaces (i.e. their paths) may 
then be obtained from their equations of variation; even when the 
group reduces to the identity, or cannot be solved for explicitly, it is 
always possible to find a basis for the differential invariants of the 
abstract space, by reducing the equations for integrability (for the 
group-equations of variation) to a finite or infinite basis. 


+ A full presentation with comprehensive bibliography to that date is to be 
found in A. D. Michal, Bull. American Math. Soc. (1939), 529-63 ; attention may 
also be called to an ingenious application by the same author in Proc. National 
Acad. Sci., U.S.A., 31 (1945), 38-43. 
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ON THE SUMMABILITY OF EIGENFUNCTION 
EXPANSIONS 


By E. C. TITCHMARSH (Ozford) 
[Received 27 May 1950] 


1. I find that, in writing my book Ligenfunction Expansions associated 
with Second-order Differential Equations, I have overlooked an impor- 
tant paper by Hilb.* The object of his paper is exactly the same as 
that of my book, viz. to prove the eigenfunction expansion of an 
arbitrary function by means of contour integration. In the case of 
series expansions, Hilb’s method would, with a little development, 
lead to the same result (Theorem 2.7(i)) as is given in my book. It 
may be claimed for my method that it can be extended to deal with 
the corresponding problem for partial differential equations,} whereas 
no such extension of Hilb’s method easily suggests itself. But Hilb’s 
method has other advantages, and I use the main ideas of it here to 
prove the following theorem. 


THEOREM. Suppose that, in the problem of the differential equation 
d? 


with the boundary condition 

$(0)cos a+¢'(0)sin a = 0, 
there is a discrete spectrum, consisting of eigenvalues Ay < A, < ...,A, > ®, 
with corresponding eigenfunctions ys, (x). Let f(x) be L?(0,00), and let 


Cn = dex. 


Then lim (x) = f(x) (1.2) 
A> x A 
wherever f(x) is continuous; and, more generally, wherever 
7 
x(n) = dt = (1.3) 


as » > 0, and so almost everywhere. 


* Math. Annalen, 76 (1915), 333-9. 
+ Proc. London Math. Soc. (3) 1 (1951), 1-27. 


Quart. J. Math. Oxford (2), 2 (1951), 258-68 
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The factor 1—A,,/A can, with a little extra analysis, be replaced by 
(1—A,,/A)®, where 8 is any positive number. Similar theorems can also 
be obtained in cases where the spectrum is not discrete. 


2. We shall suppose that sina + 0; slight modifications have to be 
made if sina = 0. Let 6(x,A), 6(x,A) be the solutions of (1.1) such that 
¢(0,A) = sina, ¢'(0,A) = —cosa, 

6(0,A) = cosa, 0’(0,A) = sina, 
where dashes denote differentiations with respect to x. Then 
W (4, 0) = (a, A)O’(x, A) —O(x, A) = 1 
for all values of x. 
Let A) = O(a, A)-+m(A)d(a, A) 

denote the solution of (1.1) which, for I(A) + 0, is L?(0,00). According 
to Weyl’s theory, this solution always exists; it is unique in the limit- 
point case, while in the limit-circle case m(A) is any point on the limit- 


circle. 
Let f(€) be any function of Z?(0,00), and let 


(x,d) = | HELE) [ HE de. 


As in Chapter II of my book, this function has simple poles at the 
eigenvalues A,. Hence, if A is real and positive, and C is a contour 
passing through A and including all the poles to the left of it, 


= | (2, (1-4) (1 (2.1) 


An<A 


It is now a question of the behaviour of (x, A) for large complex values 


of A. 
Let « = o+ir denote the value of vA such that 0 < arge < 4a 


when 0 < argA <7. Then,* for 0 < argA <7, 
A) = sina cos o(=) (2.2) 
uniformly over any finite z-interval. 
3. Consider next #(2,A) for a fixed 2. Let 6 > x be also fixed, and 
let 1,(A) be such that 
(x, A) = O(a, A)+1,(A) A) 


* Formula (1.7.3) in my book. 


t 

0 x 

= 

a 
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vanishes at x = 6. Now 


$ (a, AW (x, (w, Ap (a, A) = 8) = 1. 
Putting x = 6, it follows that 
A) = 1/(b, A). 
Hence = (3.1) 
where y(z, A) is the solution of (1.1) such that 
x(b, A) = 0, x'(b,A) = 1. 
Hence* x(x, A) = (3.2) 
« |x|? 
It follows from (2.2), (3.1), and (3.2) that, if 7 > 1, 
= sin a+ O(e-"7) + O(e"7/ 


eixe e—(2b—2)r 
t | o(—). (3.3) 
If —1 <7 <1, (2.2) gives 
$(b,A) = (3.4) 
Consider the square 
(n+1)z 


We may clearly suppose that b > 47. On the boundary of the square, 
|jcosxb| > A > 0. 

Hence, by Rouché’s theorem, if |«| is large enough, ¢(b, A) has the same 

number of zeros in the square as cos xb, viz. one. If the zero is at x = ko, 

then («—x )/¢(b, A) is regular in the square and bounded on the boundary, 

and so is bounded throughout the interior. Thus 


K— Ko 


$(b, A) 
where C is independent of x. Also x, must be real (or there would be 
another zero at the conjugate point). Hence 


1 
$(b,A)| ~ ~ 


for |r| < 4m/b, and so also for |r| < 1, by (3.4), if |«| is large enough. 
From (3.2) and (3.5) it follows that 


Pp(w,A) = (O< 7 <1). (3.6) 


* Formula (1.7.5) in my book. 


<c 


(3.5) 
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If r, denotes the radius of the circle associated with the point x == b = 
in Weyl’s theory, 


A)| = fm) )| 
< = f ae] 
0 


where v = I(A) = 207. Now 
b 


j dg = dé, 


and 
| |\cos K€ |? dé = (cos 20€+ cosh 27€) dé 
0 
1/sin2ob sinh 27b 
3\ 20 2r )> 
since sinh > 1, sin y > 
x y 7 
Hence 
>C>0 (0<7r<1) 
we dé > om (3.7) 
Hence 
=| (0<7r< 1), 
We, 
) > 1). 
Hence finally 
(0<7<1l) 
+0 )+o(—) (7 >1) 
ix sin a o 


Similar results hold uniformly if x varies over a given finite interval. 


4. We next require a result applicable to the second term in ®(z, A), 
in which (é,A) occurs with é large. As in § 2.5 of my book, 


b 
0 


U 
i 
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Taking A = p+iv, A’ = A = p—iv, this gives 


b 
Biv A) |? = A), A)} | 
0 
= W,{(x,A)+-m(A) p(x, A), O(a, A)-+-m(A) 
= W,(8, 8) +-m(A)W,(8, + |m(A) 
Let c, denote the centre of the circle C,. Then 


2 
W,(9, ¢) W,(0, 4) 


b 


= 2iv [ de {|m(A)—cy 
0 


bo 

>. 

o— > 


Hence 


b 
| Ie, der LA) f If (or, A)|2 


0 


> |? de. 
Using (2.5.2) of my book, this gives ; 
b b = 
b 0 0 


Hence, by (3.7), 


| |? dae = 
b 


—2br 
(r > 1) 
5. Consider the case where f(£) is continuous at € = x, so that, 


given e, If(E)—f(z)| <8). 
Let 


dé-+ 


= 0,+0,4+0,+0,. 


| 
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Let | O(x, ai A=14+h+h+l, 
A 


corresponding to ®,,...,0,. Let A = K? (K > 0), and take the contour 
in the A-plane so that it is symmetrical about the real axis, and so that 
its upper half corresponds to the straight lines joining K, K+iK, ik 
in the x-plane. Then 


1 
= — 
I | O(x,« ( d« 


K 
| O(z, 2+ +i) dr— 
0 


By (3.8) and (2.2), if0<7< 1, 


0 
and, if 7 > 1, 


= feve ag 
0 


(—) +40 


| 


= O(e-*/a). 


Also ®, is regular for —1 <o <1 and 7 large enough. Hence, by 
Lemma 2.11 of my book, and the last inequality, 


®, (x, x2) = O(e-*7) 


for 0 <o < 1 and 7 large enough. 
It follows that 


1 K 
0 


= O(K-1)+ O(e** K2) = 0(1) 


for a fixed § and K > oo. 


, 
| 
| 
| 
4 
| 
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Now consider ®,. We have 


x+8 


(0<7< 1), 


(7 > 1), 


by (2.2) and (4.2). By an argument similar to the above we then obtain 
I, = 0(1) as K >. 


6. Next consider ®,. We obtain, as for ®,, 
0, = (0<7<)). (6.1) 
OT 


For 7 > 1, (2.2) and (3.8) give 
= sin a+ x 
K| 


eirg e—(2b-£)7 
+8 1 x+8 
| tre | + 
(6.2) 


If —1 <o <1, this gives 


%,= 


and, using Lemma 2.11 as before, it follows that actually 
®,=O(1) <1). (6.3) 
We can also insert the term 
| ag 


on the right-hand side of (6.1) or (6.3) without altering the O-terms, 
since this term is clearly O(|«|~*). It follows that the upper half of the 
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contour to I, 


ik 

1 

K 


—(2b-22-8)r 
é T 


x 


K da 


e 


2rK g~(2b-2r-8)K 

+8 


The O-terms give of ag}+0(1) 


The main term is 


ie-KE-2) je-Ké-2) gi K(€-2) 


Since the lower half of the contour contributes the conjugate of this, 
we obtain altogether 


If we replace f(&) by f(x), we obtain 


2f (x) K(é—a) _cos | d 


(hy 
+9 | ar) 
0 
K 
i 
+0| | 
0 
— — 
™ 
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as K + 0, since 


68 


— 26? 26 
If |f(é)— = | < e, the difference made in replacing f(£) by f(x) is 


sinK(é—x) cos K(€—2) 


K(E—x)? 


K8 
= O(e) 


Hence I, > 4f(x) as K +00. It may be proved similarly that I, > } f(z) 
and this proves the theorem in the case where f(x) is continuous. 


dé 


sin 6 cos 
68 


dé = O(c). 


7. In the more general case where (1.3) holds, the last part of the 
above argument must be replaced by the following. We have 


0 <=}, 
sinK(é—x) cos K(é—z) = ( x) 


Hence 


sinK(f—a) cos K(€—z) 


@+1/K «+8 | 


= Of Kx( fal 


1/K 
The first term is 0(1), by (1.3), and 


x 


d 
= 0(K)+0(1)+0(K). 


Hence the whole expression is 0(1), and the result again follows. 
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8. It will now be shown that the above theorem includes that of 
Chapter X of my book, in which it was assumed that q(x) tends steadily 
to infinity. Let 


N 
sy = Cn Cn = S 
n=0 
Then the above theorem shows that, under the conditions stated, 
oy > f(x). Now 


N 


n=Not+1 


= 0(1), 


Cn Yn (2) 


by choosing first N, and then Ay, since > |c,,%,(x)/A,| is convergent. 
Hence 


It is easily verified that 


N-1 


Hence 


On41Ans1—FnAn 
a=? 


On An(Ansa -1) 
A? (1+A,-4/A)(1 +A,/A)(1 +An41/A) 


the rearrangement being justified since o,,/A,_, > 0. Let joy—f(x)| <e 
for VN > N,. Then this is 


A (Ay + —A, ) 
1 1 
41/A) 


A, (A, +1—An-1) 


> 


n=No+1 


| 
| 
| 
| 
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Since, for a fixed Nj, 


An(An+1—An—1) 


l< An—An-1 1 
n n O Post 
AK 2 


n=Not+1 
2 1+A,,/A A 
1 l 
~ 1+0(;), 


it follows that S > f(z). 


SOME ALGEBRAIC PROPERTIES OF CROSSED 
MODULES 


By ANNE P. COBBE (Ozford) 
[Received 31 May 1950] 


1. Introduction 

Tue term crossed module has been used by J. H. C. Whiteheadyt to 
describe groups having the algebraic properties of 2-dimensional rela- 
tive homotopy groups. In this paper some algebraic properties of 
crossed modules are derived, using the ideas of the cohomology theory 
of abstract groups developed by Eilenberg and Maclane.{ 

If ® is a multiplicative group, a crossed (8, ®)-module is defined as 
an additive, not necessarily abelian, group X, admitting ® as a group 
of operators (i.e. for which there is a homomorphism 7: ® > A(X), the 
group of automorphisms of X), together with a homomorphism 8: X > ® 
satisfying the conditions 


(Bx)y = for every x, ye X, (1.1 a) 


B(dx) = o(Bx)d-! for every xe X, GED, (1.1 b) 
where dx € X is the element obtained by operating on x with ¢. I shall 
denote this relationship by writing Bx X > ®. 

If X = BX c@is the image of X under f, it follows immediately from 
the definition that X is an invariant subgroup of ®, and we may there- 
fore form the factor group Q = ®/BX. It also follows that the kernel G 
of the homomorphism f is contained in the centre Z of X. The relation- 
ship of these groups may be displayed as an exact homomorphism 
sequence§ 

(1.2) 
where A is the identity and p is the natural homomorphism of ® onto Q. 

I shall call G the kernel of the crossed module, and distinguish between 
two cases, calling X a general or central crossed module according to 
whether G may be a proper subgroup of Z or coincides with Z. The case 


+ J.H. C. Whitehead, ‘Combinatorial homotopy (II)’, Bull. of American Math. 
Soc. 55 (1949), 453-96. I refer to this as C.H. (II). 

t Eilenberg and Maclane, ‘Cohomology theory of abstract groups (I), (II)’, 
Annals of Math. 48 (1947), 51-78 and 326-41. These will be referred to as 
C.T.A.G. (I) and (II). 

§ i.e. a sequence of groups related by homomorphisnis such that the kernel of 
each homomorphism is the image of the one preceding it. 


Quart. J. Math. Oxford (2), 2 (1951), 269-85. 
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of central crossed modules in which ® is a free group has been discussed 
by S. Maclane in connexion with operator homomorphisms of exact 
sequences.t 

The kernel @ is important in the geometrical application.t If S, T 
are arcwise-connected spaces, 7'c S, and X = II,(S,7), ® = II,(T), 
then X is a (8,®)-crossed module where f is the boundary homo- 
morphism. In this case G = where i: II,(7’) > I,(S) is 
the injection. Algebraically, G is of significance in connexion with those 
elements of ® which operate identically on X. Lemma 6 establishes this 
connexion by determining the group of these operators in terms of Q 
and BZ. 

I recall the definition of a Q-kernel§ as a pair (X, 0), where X is a group 
and @ a homomorphism of Q into the group A(X)/ J(X) of automorphism 
classes of X. It is known|| that with every Q-kernel (X, @) is associated 
a 3-dimensional cohomology class {k} ¢ H3(Q, Z), the invariant of (X, 6), 
The condition that (X,0) be the kernel of an extension} of Q is that 
{k} = 0; in this case the set of inequivalent extensions can be put into 
one-one correspondence with the elements of H?(Q,Z). These results 
are used throughout this paper. 

In this connexion we shall need to know the centre of BX. If this be 
denoted by Z*, then BZ c Z*. For, if z € Z, it follows from (1.1 a) that 
Bz is the identical operator on X. Hence, from (1.1b), with ¢ = fz, 


Bx = B(px) = $(Bx)p-. 
In general Z* is a larger group than BZ, which I shall denote by Z. 
Maclane observes{{ that every central crossed module determines a 
Q-kernel and conversely, and that this correspondence is unique up to 
equivalence. Let 8 x X > ® be a general crossed module, and let 


h: A(X) > A(X)/I(X) 


be the natural homomorphism. Using the homomorphisms p, 7 defined 
above, we obtain a Q-kernel (X, 6), where 6 is defined by the equation 


6p = hr. (1.3) 


+ S. Maclane, ‘Cohomology theory of abstract groups (III)’ (C.T.A.G. ITI), 
Annals of Math. 50 (1949), 736-61. 

+ 

§ C.T.A.G. (II) 327. 

|| Ibid. Theorems 7.1. and 8.1. 

tt (®,p) is an extension of Q by (K,@) if p: ® — Q is onto and p~(1) = K. 
Also if a € A(K) is defined by ak = dkd-, then a € O(pd). [C.T.A.G. (II) 327.] 

(TH). 
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I shall say that any Q,-kernel (X,6,) for which there is a homo- 
morphism p, of ® onto Q, satisfying (1.3) is related to Bx X >. In § 2 
I show that this relation gives rise to a unique class of Q-kernels. This 
paper is concerned with general crossed modules related to a given class 
of Q-kernels, and the two-way correspondence observed by Maclane 
appears as a corollary of the main theorems. 

We note here two extreme cases. 

(1) If no element of ® operates identically on X, then 

® = A,(X) c A(X). 
It follows from (1.1a) that B-1(1) = Z. In this case 8 x X > @® may be 
embedded in the central crossed module B x X > A(X). 

(2) If B-4(1) = 0, we may identify X and BX. Then X cQ, and © is 
an extension of Q = ®/BX by X. Any group £ of which X is an 
invariant subgroup operates on X according to (1.1) if 8 is the identity. 


2. Outline of results and definitions 

We have seen that from any crossed module 8x X >@® we can 
obtain a related Q-kernel (X,6). In Theorem 1 I obtain a necessary 
and sufficient condition, in terms of the invariant {k}, that a Q-kernel 
(X,6) be related to a crossed module B x X > ® with given kernel G. 

In the second part of the theorem I obtain a Q-kernel (X, 9), induced 
by 8 from the given Q-kernel (X,0), where 8: X > X is the natural 
homomorphism with £-1(1) = G. Then (, p) is an extension of Q by 
(X, 6); however, not all such extensions can be made to operate on X 
in the required way. I call extensions which give rise to crossed modules 
admissible and in Theorem 2 obtain the number of equivalence classes 
of admissible extensions by means of the correspondence of Eilenberg 
and Maclane between equivalence classes of extensions and two- 
dimensional cohomology classes over the centre.t{ 

Our second objective is the enumeration of equivalence classes of 
crossed modules related to a given class of Q-kernels satisfying the 
conditions of Theorem 1. Here, equivalence is defined by means of 
strict isomorphism.§ 

Crossed modules 8 x X >@® and f, x X >Q,, will be called strictly 
isomorphic if there is an isomorphism A: ® = ,, satisfying 


71,A=7 and £B, (2.1) 


t (®, p) and (®,, p,) are equivalent extensions of Q by (K,6) if A: ® = ®, is an 
isomorphism satisfying p,A = p and Ak = k for ke K. (C.T.A.G. II. 337.) 

t Ibid. (II), Theorem 11.1. 

§ Cf. operator isomorphism. C.H. (II). 
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Since B x X + © is not related to a unique Q-kernel, it is necessary to 
define also an equivalence between Q-kernels. The condition (1.3) then 
relates a strict isomorphism class of crossed modules to a unique class 
of Q-kernels. 

We say that the kernels (X,6,Q) and (X,6,,Q,) are in the same 
6-class of Q-kernels if there is an isomorphism 7: Q ~ Q, satisfying 


= 0. (2.2) 


Then Q-kernels related to strictly isomorphic crossed modules are in the 
same 6-class. 

For let 7: Q = Q, be defined by np = p,A, where A is given by (2.1). 
Then, since p, p,; are homomorphisms onto, and A is an isomorphism 
onto, 7 is a homomorphism onto. Let ng = 1 for some gE Q. Then 


1 = = = pi Ad 


for some ¢€®, and so Ade f, X. From (2.1) it follows that ¢ € BX; 
therefore ¢ = pd? = 1, and 7 is an isomorphism. 
From (1.3) and (2.1) it follows that 


0, np = 9, = hr, A = hr = Op. 


Hence (2.2) is satisfied. 

We need also a definition of equivalence of extensions of kernels 
(K,@,Q) and (K,6,, Q,) in the same 6-class. The extensions (®, p) of Q 
by (K, @) and (®,, p,) of Q, by (K, @,) will be called 6-isomorphic if, given 
7: Q = Q, where 6, » = 9, there is an isomorphism A: ® ~ Q, satisfying 


pPiA=np and Ak=k for every ke K. (2.3) 


The case Q = Q, and n = 1 reduces this to the equivalence defined by 
Eilenberg and Maclane. 

For simplicity I prove two lemmas showing that no generality is lost 
by taking 8 and X as fixed (Lemma 3) and Q and 6 as fixed (Lemma 4). 
The problem then reduces to finding strict isomorphism classes of crossed 
modules from the admissible extensions of a fixed group Q by a fixed 
Q-kernel (X, 

It follows from the above working that, if Bx X > ® and Bx X > 9, 
are strictly isomorphic and related to the same Q-kernel (X,6), then 
(®, p) and (®,, p,) are 6-isomorphic extensions of Q by (X,9); 7 is now 
an automorphism of Q satisfying 07 = 0. This is a stronger condition 
than that of (2.3), which it implies; (2.3) demands only that 6 = 0. 
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In §9 I find the condition that two extensions be 6-isomorphic, 
and restate Theorem 2 in terms of isomorphism classes of extensions. 

Finally, to complete the enumeration, I find the number of strict 
isomorphism classes of crossed modules obtainable from a fixed iso- 
morphism class of admissible extensions of Q by (X, 6). 

In order to simplify the statement of the theorems, I show in the 
next section how the homomorphism f can be used to obtain two 
relative cohomology sequences in terms of which the results may be 
stated. 

For reference I give a list of the symbols used throughout the paper: 

A(X) is the group of automorphisms of X; 

J(X) is the group of inner automorphisms; 

C{x] is the inner automorphism of X induced by 2; 

h: A(X) > A(X)/ I(X) is the natural homomorphism; 

® operates on X through the homomorphism 7: ® + A(X); 

B: X + © is the crossed module homomorphism and 7(fx) = C[2]; 

Q = 0/BX and p: © > Q is the natural homomorphism; 

6: Q + A(X)/1(X) is defined by 6p = hr; 

X = BX, Z = BZ, where Z is the centre of X; 

Z* is the centre of X and Z* = B-(Z*). 


3. Relative cohomology sequences 
Let T, I’ be two additive abelian groups admitting a group Q of 
operators, and let 


be a homomorphism with kernel I. If C"(Q,T) is the group of n-co- 
chains on Q to T, and fe C"(Q,T), then A induces a homomorphism 


r’: T) + CX(Q, 

Let A be an operator homomorphism, i.e. A(qa) = q(Aa) for ae T. 
Then A’8 = 8A’, where 5 is the co-boundary operator; i.e. when 
fe crQ,T), (8f) C"+(Q,T) is defined by 

n 


+ (— 


If i’: 0n(Q, CX(Q,T) 
is induced by the identity map i: Tj > I, the kernel of 2’ is i’C"(Q, T)). 


3695 .2.2 T 


1 
: 
= 
= 
hi 
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For simplicity I shall take 7’ as the identity and write \’-1(0) = C"(Q,J}). 
Since [, is closed under the operations of Q, therefore 


We have, therefore, a system of co-chain groupst 
= {C"(Q,T)} 


with C(T) c C(t) and A’C(T) = c C(I’). We may define the 
groups and [}) of relative co-cycles and co-boundaries 


and dfe c+ 
fe if and f—she C(Q,T) for some 
he C"-(Q,T). 
The relative cohomology group is then defined by 
= 2“ 
With these definitions A induces 
A: H"(T) > AAP), 
jy: & HD), 
where A is a homomorphism onto and j) is an isomorphism. 
We obtain the two exact homomorphism sequences 
—> H(T,) Hr) A> —, (3.1) 
—> Hr) —ar) Bear) —> (3.2) 
where i, i* are injections and j* is induced by the natural homomorphism 


I’ + 
In our application we take fT = Z, Ty = G, I’ = Z*, andA = AZ. 


4. THEOREM 1. (A) Let Bx X >® bea crossed module with kernel G. 
If {k} is the invariant of a related Q-kernel (X,0), then 
{k} € iH(Q, @). 


(B) Conversely, let (X, 6) be a Q-kernel with invariant {k}, and let G be 
a subgroup of the centre of X such that 0(q)G = G for everyqe Q. If 


{k} @), 


there is a crossed module Bx X >@® with kernel G, related to the given 
Q-kernel (X, @). 


+ Kelley and Pitcher, ‘Exact homomorphism sequences in homology theory’, 
Ann. of Math. 48 (1947), 682-709. 
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To prove the first part of the theorem we consider the exact sequence 

The conditions (1.1) are those imposed by S. Maclane on the exact 
sequences with which he associates a 3-dimensional cohomology class 
{y} ¢ H°(Q, @), called the obstruction class} of the sequence. Theorem 1 A 
follows immediately from the following 


Lemma 1. Let G6>X 404% Q be an exact sequence satisfying (1.1), 
with obstruction class {x}. The related Q-kernel (X,0) has invariant {k}, 
where {k} i{x}. 

Let 7, h, 6 be the homomorphisms defined above. We first construct 
the obstruction x(¢,, 42,93) defined on Q with values in G. 

For every q € Q, choose a representativet ¢(q) € ® such that 


= 


These representatives determine a factor set 2%(q,,q,)¢BX = p-(1) 
satisfying the associative condition 


$(91)% (Go, Ys){$(91)} 929s) = (G1 Ys)- 
For every £(q;, 92) choose x(q,,q2) € X such that 
= 
Using (1.1b), we may write the above associative condition 


The two sides of this equation determine an element 
42.93) € G = B-*(1), 
satisfying 


Different choices of representatives determine co-cycles cohomologous 
to x, and the sequence determines a class {x}. 

To find a particular k for the Q-kernel (X, 6) we must « xoose, for each 
q € Q, an automorphism a(q) € 0(q). If ¢(q) is the element chosen above, 


we may choose a(q) = 76(q) 


since ha(q) = hrd(q) = = 9(q). 


C.T.A.G, (IIT). 
t In choosing representatives, e.g. $(q) € ® for gE Q we may take ¢(1) = 1. 
This choice is assumed throughout the paper. 
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If we apply 7 to the equation 


we obtain = Cl x(q, G2) (4.2) 
where x(q, 42) is the element chosen above. 
The equation defining & is then 
Is) 9293) = + +2(91 929s). (4.3) 


Comparison with (4.1) shows that, with this choice of representatives, 


929s) = Ys) 
whence {k} = ify}. 


5. Before proving the converse, I prove a general lemma on Q-kernels. 


Lemma 2. Let (X,6) be a Q-kernel with invariant {k}, and let u: X + Y 
be a homomorphism onto a group Y. 

If, for every «€h-*0Q c A(X), af{u-1(1)} = p(1), then p induces a 
Q-kernel (Y,0) with invariant {uk}. 

Let A,(X) = h-10Q c A(X). For each «€ A,(X), define wa = & by 
Op = pa. 

Since aa € p~1(1) whenever x € .~1(1), it follows that & is single-valued 
and hence an endomorphism of Y. 

Since, also, ax € w~*(1) implies x € p-1(1) for all « € A,(X), therefore 
ae A(Y). 

For a4, % € Ap(X) and xe X, 


= p(x, 2) 


= & 

= 
and, since &, &, %a@ are automorphisms this implies = dp. 
Therefore w: A(X) > A(Y) 


is a homomorphism. Moreover, wC[x] = C[u«], and, since yp is onto, 
therefore wI(X) = I(Y). Accordingly w induces a homomorphism 
@: 0(Q) > A(Y)/I(Y), where ah = hw with h: A(Y) > A(Y)/I(Y) the 
natural homomorphism. 

By defining 6 = 0, we obtain a Q-kernel (Y,6), where ha = 6q 
implies 


= hwa = = = 6q. 
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By applying w to (4.2) and p to (4.3), we see that the invariant of 
(Y, ) is {uk}. 


6. To prove the second part of Theorem 1, define 
p:X>X = X-—G 
to be the natural homomorphism. 
The conditions of Lemma 2, with f for » and X for Y, are satisfied, 
and we obtain an induced Q-kernel (X, 9) with invariant {k} = {Bk}. 
If {k} e iH3(Q, @), we may choose k(q,,q2,93) to have values in G. 


We may then chooset a(q) € 0(q) and x(q,,q.)€ X for each ge Q, so 
that (4.2) and (4.3) are satisfied. 


Since k(q,, 42,3) € G for all q € Q, therefore k = Bk = 0 for this choice 
of k. Hence {k} = 0. This is the condition that the induced kernel 
(X, 8) be extendible. 
Let (®, p) be the extension of Q by (X,6) whose elements are pairs . 
(@,q) with p(%,q) = q, and multiplication defined by 
2) = (€ ) Fo} qz)- 
By choosing 4(q) = (1,q) and identifying (#,1) with ¢ ¢ X we have 
= V2) (91 
= a(q)e. 
We make ® operate on X, defining 7: ® > A(X) by 
= Clx]a(q). 
Then 7 is a homomorphism. For 
= Cla, 92) 2) 
= 
= Cla 
= 7k, $(91)}7 {Xe $(2)}- 


The operations of ® on X satisfy (1.1), and we have a crossed module 
Bx X > © related to the original Q-kernel (X, @). 


7. Admissible extensions of (X, 6) 
In order that 8 x X > © be related to the Q-kernel (X, @), (®, p) must 
be an extension of the Q-kernel (X, 6) induced by £. In the last section 


+ C.T.A.G. (IT), Lemmas 7.2 and 7.3. 
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one such extension was obtained. We seek the condition that another 
non-equivalent extension be admissible. 

Let a(q) € 0(q) be chosen for each ge Q. We shall assume that this 
choice of a(q) remains fixed. Let (®,,p,) be the extension already 
determined. Choose ¢,(q) € pj 4(1) for each q so that 


71 = 
and let #,(¢,,q.) be the factor set corresponding to ¢,(q). Then «(q) 
isfies (4.2) with Fa 
92) = Y2)}- 
Let (®,, p.) be any other extension, and ¢,(q) € pz 4(1) be chosenf so 


If #,(9,, 92) is the corresponding factor set in ®,, then 
92){F1 (G1 = Ge), (7.2) 


where 2z*(q,,92) is a co-cycle in C*(Q,Z*). The two extensions are 
equivalent if and only if this co-cycle is a co-boundary, and it is by 
means of the cohomology class {z*} « H?(Q, Z*) that the correspondence 
between equivalence classes of extensions and elements of H?(Q, Z*) 
is established. 

Suppose that (®,,p,) is an admissible extension and the resulting 
Bx X +9, is related to the Q-kernel (X,@). There is a homomorphism 
T,: ®, > A(X) satisfying (1.1) and (1.3). Therefore 

hrs $x(0) = = = hag). 


We may then assume that ¢,(q) is chosen so that 
T2$2(7) = = 719, (9). (7.3) 


This choice of ¢,(q) evidently satisfies (7.1). 
We now have 


C[%2(91, 92) ] = 
= 
= (91 92) 
= V2) 
= %2)} 
= 
Hence € Z, 


¢ C.T.A.G. (IT) 337, 338. 
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and from (7.2) we have that, for this choice of representatives, 


2* (41,42) € Z 
for all g,, € Q. This implies that 
i*H*(Q, Z). (7.4) 

Conversely, let (®,,p.) be an extension for which the co-cycle z* 
defined by (7.1) and (7.2) satisfies (7.4). By a suitable choice of 
$(q) € pz (q) we may take z*(q,,q2) to have values in Z. 

Define 7,: ®, > A(X) by 7,{%4,(q)} = Cla]a(q). Since 

C[ = G2) ] = 92) 

7, is a homomorphism satisfying (1.1), and the extension is admissible. 

We thus obtain 

THEOREM 2. The equivalence classes of admissible extensions of Q by 


the induced Q-kernel (X , 8) are those corresponding to elements of H®(Q, Z*) 
which lie in itH%(Q, Z). 


8. The following lemmas show that there is no restriction imposed 
by considering only one induced Q-kernel (X, 6). 

Lemma 3. Let X ~ X—G be a given homomorphic image of X, and 
let By: X -> X be a given homomorphism onto X with By (1) = G. Inany 
class of strictly isomorphic crossed modules B x X > ® with given X and 
G, there is always one such that BX = X and B = By. 

Let B,x X +, be any crossed module of the given class, and let 
A: B, X > X be defined by NB, x) = Bow. 


Then A is an isomorphism onto X. 
Let ® be the aggregate of elements of X together with those ¢ €®, 
which do not lie in 8, X. A one-one map A: ®, > ® is defined by 


Ad = ¢ for ¢ not in B, X, 


= Boe. 
If multiplication in ® is defined by 
(Adi = A(P $2), 


then © is a group, and A an isomorphism onto ®. 
We define 7: ® > A(X) by 


TA = 17; 
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and obtain (Box) = = = 
(Ab )(Box)(AP)* = 
= NB, (¢x)} 
= Bol(Ad)a}. 
The conditions (1.1) are therefore satisfied by B, and ®, and 
ByxX >® 


is in the given class and satisfies the requirements of the lemma. 

We wish to show also that, in constructing all the strict isomorphism 
classes of crossed modules related to a given 6-class of Q-kernels, there 
is no loss in considering only one Q-kernel of the class. 

In § 2 it was shown that, if (X,0,,Q,) and (X, 6,, Q.) were related to 
strictly isomorphic crossed modules, they were in the same 6-class. 
Conversely I prove 


Lemma 4. Let Bx X +® be related to (X,6,,Q,) and let (X,6,Q) be 
a Q-kernel in the same 6-class. Then Bx X >® is related to (X,6,Q). 

Since Bx X > © is related to (X,0,,Q,), there is 2 homomorphism 
p;: D > Q,, onto Q,, such that 6, p, = hr. 

Since (X,0,,Q,) and (X,6,Q) are in the same 6-class, there is an 
isomorphism 7: Q, ~ Q such that 6 = 4). 

Define p: ® > Q by p = np). Since p, is onto and 7 is an isomorphism, 
therefore p is a homomorphism onto Q. Also 

Op = Onp, = 9, py = hr. 

Consequently (1.3) is satisfied and (X, 6, Q) is related to Bx X >. 


9. In § 2 it was shown that, if 8x X > ® and Bx X >, are strictly 
isomorphic, then (®, p) and (,, p,) are @-isomorphic extensions of Q by 
(X,@). This implies that the extensions are 9-isomorphic. For, if 
n € A(Q) satisfies 6y = 0, then, in the notation of § 5, 

= = = 8. 
The converse is not necessarily true. 

It is clear that extensions (H,, p,) and p,) of Q by (K,@) which 
are equivalent are 6-isomorphic. I seek a condition that non-equivalent 
extensions be 6-isomorphic. 

Let [ be an abelian group on which Q operates, and let 7 € A(Q) be 
such that (nq)c = ge for all ce T, ge Q. Define an automorphism 7’ 
of C"(Q,T) by 

Since = 75, induces an automorphism 7* of H"(Q,T). 
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We obtain an equivalence relation between elements of H"(Q,T) by 
setting {c,} = {c} whenever {c,} = *{c}, where y* is induced by some 
n € A(Q) as above. 

If [' is the centre of K, and » is a 0-automorphism of Q, i.e. 6n = 8, 
then (nq)c = ge for allce T,qgeQ. 

We denote by H3(Q,T) the set of equivalence classes of elements of 
H?(Q,T) under the above equivalence relation induced by 6-auto- 
morphisms of Q. We then obtain a condition for isomorphism of 
extensions. 


Lemma 5. The 6-isomorphism classes of extensions of Q by (K,6) are 
in one-one correspondence} with the elements of H}(Q,T). 


We assume, as before, that «(q) € 0(q) is chosen for each g € Q, and 
remains fixed. If, therefore, »¢A(Q) is a 9-automorphism, since 
4(nq) = 6(q), then a(nq) = a(g). Let (£;, p;) (¢ = 1, 2) be 6-isomorphic 
extensions of Q by (K,@). There is then an isomorphism A: FE, = E, 
satisfying p»A = np, for some 7» € A(Q) with Oy = 0. 

Corresponding to the choice of a(q) € 6(q) we may selectt elements 
42) € K so that the resulting 3-co-cycle k(q,, 42,93) = 1. Let 
be the group of pairs (k,q) with multiplication defined by 


o,(k,q) = q and f,(q) = (1,q). Then f,(q;,q2) is the factor set. 
The extension (Z,, p,) may be described§ by choosing e,(q) € pz 1(q) to 


The corresponding factor-set for is then Where 
€ the centre of K. 
Let a function f,(q,,¢2) be defined for all q,, ¢.€ Q by 


192) = 12): (9.3) 
and let (F,,c,) be defined in the same way as (F,,0,) with f,(q,,q2) in 
place of 2). 

We can describe (E,,p,) by choosing e,(q) € pz to satisfy an 
equation similar to (9.2) and having the function ¢2(q,,42)fo(41,%2) as 
factor set. 


+ This correspondence is analogous to that between equivalence classes and 
elements of H2(Q, Z). 

C.T.A.G. (II) 333. 

§ Ibid. (II), Lemma 11.2. 
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Since A: EZ, ~ E, satisfies p.A = np,, there are elements k(nq) € K such 


h 
ey(g) = (9.4) 
Since for this y € A(Q), (nq) = a(q) and Ak = k for every k € K, we 
have from (9.2) and a similar equation 
= a(nq)k 
= a(q)k 
= 
= 7, 
whence, using (9.4), k(nq) € T for all q. 
Applying A to the equation 
we have 
from which it follows that (writing [' additively) 
92) = 192) + 092). (9.5) 
Hence {c.} = n*{c,}. 
Conversely, if this relation holds for some 0-automorphism 7, there 
are elements k(nq) ¢ T, with k(1) = 1, satisfying (9.5). 
Let A: E, > FE, be defined by 
= for all e K, qe Q. 
Then A is an isomorphism onto, satisfying p.A = np,, and the extensions 
(Z,, p,) and are 0-isomorphic. 
We may now rewrite Theorem 2 in terms of isomorphism classes of 
extensions. 


THEOREM 2A. The 6-isomorphism classes of admissible extensions of Q 
by (X, 6) are in one-one correspondence with the classes of H}(Q, Z*) which 
lie in i*H}(Q, Z). 


Using the notation of the last lemma, with K = X, we have 
42) € X, Z* (i = 1, 2). 
As in § 4, the extension (F,,0,) is admissible, where 7, is defined by 
71(2,q) = C[x]a(q). 
Let x,(9,,%2) be chosen in B-1/,(q,,q2) and x2(q,,q2) be defined by 
(915 192) = 


th 


| 

q 
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Since, when a(q) € 6(¢g) are chosen automorphisms of X, if 6m = 9, 
then «(nq) = a(q) = a(7-1q), it follows that 
G2) |] = Ye) } 
= 
= 
Therefore 
The extension (F,,c,) is therefore admissible where 7,(%,q) = C[«]«(q). 
Let (Z,, p,;) be an admissible extension. From Theorem 2 this implies 
that {c,} = i*{Z,} for some {2,} e H2(Q, Z). We may assume 
92) = Y2)- 
Define 2,(91, 192) = for all q,, Q. From (9.5) we have 
Cy = 2,—8k where k C\(Q, Z*). Hence 
{co} = = = 
From Theorem 2 it follows that (£,, p,) is admissible. 


We have, therefore, that an extension which is 6-isomorphic to an 
admissible extension is itself admissible and the theorem follows. 


10. To complete the enumeration, we must determine the number of 
strict isomorphism classes of crossed modules that can be obtained from 
a given 6-class of admissible extensions. 

THEOREM 3. The classes of strictly isomorphic crossed modules derived 
from a given 6-isomorphism class of admissible extensions are in one-one 
correspondence with 1*-1(0). 

Let (®,,p,) and (®,, p,) be 6-isomorphic admissible extensions. Let 
$,(q) € pz and 4,(q) € pz be chosen so that (7.1), (7.3) are satisfied. 

Since the extensions are 9-isomorphic, there is a A: ®, = ®, such that 
pA = np, for some 6-automorphism 7» € A(Q). Then 


A$ (q) = (10.1) 


where, as for the function k(q) defined in (9.4), z*(q) € Z* for all q. 
Since 7 is a 6-automorphism, we have from (7.3) 


71 = a(¢) = = T2$2(79)- 
Hence the factor sets and ¥,(q,,q2) satisfy 
and therefore 92) = 192)%1(%1 V2) (10.2) 
where 2(q,,q2) is a co-cycle in C*(Q, Z). 


} 
| ‘ 
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It follows from the equation 


using (10.1) and (10.2), that 


2091» 192) = —(82*)(91, 092)- (10.3) 

The choice of ¢,(¢) and ¢,(q) may be altered by multiplying by a factor 
in Z, and the given conditions will still be satisfied. This would alter 
the co-cycle 2Z(¢,,q2) to one differing from it by a co-boundary in Z. 
Hence the two extensions determine an element of H2(Q, Z) which lies 
in i*-1(0). 

If the crossed modules 8 x X > ®, and Bx X >, are strictly iso- 
morphic and} A satisfies 7,A = 7,, it follows from (10.1) and § 11 below 
that z*(q) has values in Z. Hence, from (10.3), 2(q;, 4) is a co-boundary 
in Z. 


Conversely, let = (8£)(41; Ye) 


( 
for some 1-co-chain £ ¢ C1(Q, Z). Choose new representatives 


$2(9) = € 
with the factor set = {(8Z)(q,, 
From (10.2) it follows that 


192) = G2) 
and 72 = = 71 $1(9)- 
Let A’: ©, > ®, be defined by 
= $2(09), NE =F. 
Then A’ is an isomorphism satisfying 7,’ = 7,. Hence the crossed 


modules are strictly isomorphic if and only if {Z} = 0, and the theorem 
follows. 


11. Identical operators 


The algebraic significance of G the kernel of B x X + ® is connected 
with the identical operators of ®, i.e. elements such that dz = x for 
every X. These form a subgroup = 7~-1(1). 


Lemma 6. Let Q, = 0-(1)c Q. The group ®, of identical operators is 
a central extension of Q, by Z. 


{ If the strict isomorphism were in virtue of another isomorphism 2’, this 
would satisfy p,A’ = 7’p,, for some 6-automorphism 7’ € A(Q). The argument 
would be the same with 7’ in place of 7. 
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If Ba € Do, it follows from (1.1 a) that 
x+y—ax = (Bx)y = y for all ye X, 


and so 2 € Z. The intersection of ®, and X is therefore Z. 
If dy € Do, it follows from (1.1 b) that 


Bx = = 
and so X commutes with ©). Therefore Z is in the centre of ,. 

Since 6p¢, = hrdy = 1, therefore pdy € Yo, and p induces a homo- 
morphism py: > Qo. Since p-(1) = X, therefore p> 1(1) = Z. 

Also, given any qd) € Qo, there is some ¢ € ® such that pd = gy. Then 
hrd = Oph = 1, whence ro € I(X). Let 74 = C[x]. Then 

$(Bx)* €r-*(1), 
and ¢ = ¢,(Bx) for some Dp. Since pydy = pd = qo, the induced 
homomorphism p, is onto, and ®, is a central extension of Q, by Z. 

It follows that the existence of identical operators depends only on 
the two homomorphisms @ and 8. Those arising from the kernel of @ 
impose no restriction on X, since any Q-kernel (X,6) is related to a 
central crossed module. The existence of those arising from B depends 
on whether 8-1(1) is a proper subgroup of Z, and thus reduces to the 
problem of Theorem 1. 


12. Central crossed modules 

The identification of central crossed modules and Q-kernels observed 
by Maclane follows immediately from the theorems. If G = Z, the con- 
ditions of Theorem 1 are always satisfied; a Q-kernel therefore always 
gives rise to a central crossed module. In this case Z = 1 and H*(Q, Z) 
contains only one element. Theorem 2 then shows that there is only 
one admissible extension, and Theorem 3 that this extension leads to 
only one class of crossed modules. The correspondence is therefore 
unique. 

I should like to thank Professor Whitehead for the many helpful 
suggestions he has made during the preparation and writing of this paper. 
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ON MATRICES WITH POSITIVE ELEMENTS 
By HAZEL PERFECT (Swansea) 


[Received 25 July 1950] 


1. Introduction 

OsKAR PERRON (1), (2) was the first to prove that a matrix A, all of whose 
elements are real and positive, has a real and positive characteristic root 
x which is simple and greater in absolute magnitude than every other 
root; and that the elements of the adjugate matrix of y/—A are all real 
and positive: x is called the maximal root of A. The same result was 
later proved again by Frobenius (3) and extended in one or two direc- 
tions. His proof, in contrast to Perron’s, is purely algebraical and is 
more direct. A similar result for irreducible matrices with non-negative 
elements has been given by Frobenius (4) and, very recently, by 
H. Wielandt (6). 

The present paper describes a simple approximation method whereby 
the maximal root of a positive matrix can be approached. Perron was 
interested in this aspect, and in fact his proof of the existence of the 
maximal root depends on proving that 


for all i, k, where A’ = (a{?), etc. Theorem I below, in its present form, 
is due to Collatz (5), but in essence it is containec in the paper (3) by 
Frobenius. Theorem IT has not been stated before, but Perron uses a 
similar type of proof to prove a somewhat analogous inequality. The 
proof of Theorem IIT has been modelled on a method of Perron’s. The 
concluding geometrical supplement on symmetric matrices with positive 
elements is new. 

I am grateful to Dr. R. Rado for his kind encouragement and helpful 
criticisms. 


2. Proofs of Theorems I-III 

I use wu to denote a column vector whose elements are real and positive, 
but are otherwise arbitrary. A column vector with elements 2,, %»,..., 2, 
is denoted by {,, x2,..., 2,}. The particular column vectors 


{1,0,0,...,0},  {0,1,0,..., 0}, 


are denoted by ¢,, ég,.... 


Quart. J. Math. Oxford (2) 2 (1951), 286-90. 
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THEOREM I. J'he maximal root x of A is included between the greatest 
and the smallest of the numbers = (i = 1, 2,..., n). 


Proof. According to Perron’s theorem quoted above there exists a 
vector z with positive elements which is the pole corresponding to the 
characteristic root x of A’; so that A’z = xz. Write 


M = diag (py, 
Now = 2’Au = xz'u, 
i.e. z'(xI—M)u = 0, 


from which it follows that the n numbers x—p, (i = 1, 2,..., m) cannot 
be all positive or all negative. 


THEOREM II. Let r be a given positive integer, and let pu” be defined by 


Let max = = ming? = 
Then pO > (6 = 1, 2...., 2). 


where J,, J, are obtained from the unit matrix I by respectively multi- 
plying and dividing the elements of the sth row (s = 1, 2,...,n) by e, A’u. 


Then 
and (e; AZ,){1, 1,..., 1}. 
Moreover ¢, AL{I1, 1,..., 1} = 
Hence > pl) > for all i. 


THEOREM III. The sequence un” increases and the sequence ji” decreases, 
in the wide sense, to a common limit which is the maximal root of A. 


Proof. By Theorem IT, 
< porn < 


| 
| 
A 
re 
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Hence there exist limiting values 


lim = pg, lim porn fi; 
r—oo ro 
and pcp. 


Since x is known to lie between u™ and ~” (Theorem I), it remains to 
prove that » = ji. 
Now, if « denotes an arbitrarily small number, then, for sufficiently 


large r, > > 


Therefore 


e; "Attly 


and, since every element of A is real, positive, and finite, every element 
of the sth row (column) of A is less than some positive number N multi- 
plied by the corresponding element in the pth row (column). 
Ae, Aru] 


h 
Therefore Ae, Aru] 


It follows that D—p > 
This is true for all i, p; hence, in an 


Therefore 


Since this is true for all positive e, 


b 


tic 
| 
| 
to 
| 
— i 
] 
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Note. Inthe more restricted case when A is also symmetric, Theorem 

III can be proved very simply by considering an orthogonal transforma- 
tion P’AP. 


3. Geometrical supplement on symmetric matrices with positive 
elements 

The fact that the numbers yp!” (i = 1, 2,..., n) enclose the maximal 
root x of A more and more closely as r > 00 links up rather satisfactorily 
with the following geometrical result proved here for a positive and 
symmetric matrix. (Certain numerical examples would lead one to 
suppose that the result is true for any positive A.) 

THEOREM. The projection of the unit vector A"u/|A*u| on the pole, y say, 
corresponding to the maximal root x of A increases as r increases, and tends 
to the value 1 as r tends to infinity, i.e. the angle between the vectors A'u 
and y tends to zero as r > ©. 

(The projection of a vector x on a second vector z (+ 0) is defined to 
be z’a/|z|; it is well known that 

2'x 
with equality only when 2/|x| = z/|z|.) 
Proof. Consider the square of the number 
Since A is symmetric, there exists an orthogonal matrix P such that 
P'AP = diag(x, p,..., w). 
Also P’A?"P = diag( ?",..., for every positive integer r. More- 
over, since x is the maximal root of A and since it is simple, 
Put a= PU, gy = PY. 


Then, since Ay = yy, and P is orthogonal, 
APY = xPY and P’APY = 


ie. diag(x, p,...,w)Y = xY; and, since x ¥ p,..., x # w, it follows that 
=... = Y, = 0. Since Y 0, Y, is not zero. 
Further A'y = x’y, so that y’A’ =,’y’. 


3695 .2.2 U 
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Accordingly 


U3 


PUR... + U3, 
But U, £0 since y’u = Y,U, > 0, and therefore this expression in- 
creases steadily as r increases, and tends to 1 as r > oo. 
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THE APPLICATION OF VIRTUAL SOURCE 

DISTRIBUTIONS TO PROBLEMS IN THE 

LINEARIZED THEORY OF SUPERSONIC 
FLOWS 


By J. F. W. BISHOP (London) 


[Received 17 August 1950] 


Summary 

THE generalized source-sink method of the linearized supersonic flow 
theory is used to solve certain problems in two-dimensional and axially 
symmetrical nozzle flows, and in aerofoil theory; the term ‘virtual 
distribution’ is used to differentiate the method from the familiar 
equivalent-layer theory, which is seen as a special case of the method 
of virtual distributions. It has not been found possible, in general, to 
find the distribution required to solve a given problem, and the inverse 
approach, by means of which the flow patterns corresponding to a 
specified distribution are determined, is used. 

The results include the solution of a specific axially symmetrical 
nozzle flow developed in terms of elliptic integrals, the solution of a flat 
inclined plate with a hyperbolic leading edge, and the discussion of a 
family of symmetrical wings, one member of which is examined in detail. 


1. Introduction 

The use of sources, sinks, and doublets is an established feature 
of classical hydrodynamics, and the extension of the method to the 
linearized theory of compressible flows follows naturally from the formal 
similarity of the potential equations. 

Several writers have considered problems in the supersonic region 
using equivalent-layer distributions; the object of this paper is to 
examine certain flow patterns derived from the more general class of 
arbitrarily located distributions of which the equivalent layers consti- 
tute a special case. 

It will be shown that the method is applicable both to aerofoil theory 
and to channel and nozzle flows. 


2. The linearized theory 
Investigation of the irrotational, isentropic flow of an inviscid, non- 
heat-conducting perfect fluid leads to the consideration of a non-linear 
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second-order partial differential equation for the potential; if considera- 
tion is restricted to flows in which only small modifications of a uniform 
stream occur, the technique of small perturbations may be employed. 
If a steady stream of (supersonic) velocity U, sonic velocity a, and 
of density p is taken to flow in the positive x-direction, the Prandtl- 

Glauert equation for the potential is 


(2.1) 


and since the perturbation potential ¢ is related to the total potential © 
b 
f © = Ux+4, (2.2) 


the perturbation potential (and its derivatives) will satisfy the linearized 
equation. The assumption made in the derivation of (2.1) is that the 
quotients of the perturbation velocities (uw, v, w in the x, y, z directions 
respectively) by the stream velocity U must be negligible compared 
with unity. 

The differential form of Bernoulli’s equation, 


leads to the result p= —pUu (2.4) 


for the excess pressure at a point. 


3. Formulation of the problem 

An isolated point disturbance is propagated within a cone lying 
downstream from the point with semi-vertical angle cot-1/(U?/a?—1). 
This ‘Mach’ cone is therefore a surface upstream from which the flow 
is undisturbed, and within which the basic stream is subjected to small 
perturbations. All slender disturbing bodies will have a similar asso- 
ciated surface known as a Mach front, which may be found as the 
envelope of a Mach cone whose vertex maps the body. 

The two regions of the flow separated by this Mach front are con- 
sidered separately, and the solution to the problem is non-analytic. 
Some extra fact is therefore needed to determine the flow pattern 
uniquely, and the physical limitation to the propagation of (small) 
disturbances supplies the required condition; in practice this is usually 
imposed by the application of the shock-wave relations at the Mach 
front. 
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A solution of the linearized equation derived from a purely mathe- 
matical approach is therefore required to satisfy both the boundary 
conditions on the body and the shock-wave relations at the Mach 
front. 

Solutions derived from source, sink, and doublet distributions embody 
the physical limitations to the propagation of disturbances, and the 
problem is reduced to isolating a distribution which will satisfy the 
boundary conditions on the body. 


4. Channel and nozzle flows 

The two-dimensional linearized theory of channel flows is very simple, 
whether approached from the general solution to the linearized equation 
or from the linearized characteristics method. Axially symmetrical and 
general three-dimensional flow problems are more difficult, and, although 
the axially symmetrical flow problem has been examined by using the 
general solution in terms of Bessel functions both for a free jet (1) and 
for a quasi-cylindrical tube (2), the source-sink method does not appear 
to have been used. 

Virtual source distributions are applicable to both two- and three- 
dimensional problems. It has not been found possible, in general, to 
determine the density function required for a specified flow, and the 
inverse problem is considered. The two-dimensional flow is considered 
first since the results, although easily obtainable otherwise, indicate the 
general procedure to be followed. 

Associating auxiliary coordinates €, », ¢ with the x, y, z directions 
respectively, virtual distributions lying in the (»,¢)-plane will be 
examined. If the two-dimensional flow is independent of y, the source 
density will be a function of { only, and the potential at P (x,y,z), due 
to a source distribution of density f(¢) will be 


axe dnd 


where is the region of the (n, f)-plane lying within the forecone of P 
(i.e. a Mach cone lying upstream from P). 


0 <a), 
>a), 


As anexample let = 


i.e. consider a uniformly dense source plane containing a slit of breadth 
2a parallel to and symmetrically disposed about the axis. 


al 
| 4 
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The perturbation potential must be evaluated in three regions (Fig. 1), 
and it may be shown to have the values 


regionI: ¢@=0; 


— ed 7 A o¢ 

region III: ¢ = 


7~ 
I 
L 


~ 
~ 
~ 


Fic. 1. Integration regions for a two-dimensional nozzle flow. 


If the two streamlines arriving at the slit edges are taken to bound a 
channel and A = + Ua/z, this virtual distribution is seen to correspond 
to a divergence (or convergence) through a wedge section of angle «, 
with subsequent reversion to a broadened (or narrowed) uniform stream 
as shown in Fig. 2. By a suitable selection of the density function a 
channel of any shape (not necessarily symmetrical about the 2-axis) 
may be treated. 

The axially symmetrical and general three-dimensional nozzle flow 
problems can be treated in an exactly analogous manner by making an 
orifice in the plane and again considering as the initial stream a stream 
tube which does not intersect the source plane (not necessarily the tube 
of which the orifice is a section). 


Il 
j 
: 
~ 
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A discussion of the type of density function from which useful results 
may be expected will be given when an example has been considered. 


4 
—U> U+ y 


Fic. 2. Divergent and convergent two-dimensional flows 
corresponding to a uniform slit source plane. 


5. An example in axially symmetrical nozzle flow 

The flow pattern produced by a source plane is dependent both on 
the density of the source distribution and on the geometry of the orifice; 
in this section the simplest available three-dimensional case will be 
examined. Consider a density function 


i.e. consider a source plane of density »?+¢?—a? containing an orifice 
of radius a centred at the origin. 

The associated flow will be axially symmetrical and the stream lines 
for a nozzle flow which is initially cylindrical may be examined in any 
meridional plane; taking this to be the plane y = 0 and introducing a 
radial coordinate r to replace z, the perturbation potential at P (x, r) is 


ly (5.2) 


where > is the region of the source plane lying within the forecone of P 
(this is equivalent to associating rectangular coordinates y and { with 
the plane «= 0 where the ¢-axis is parallel to the radius vector 
through P). 
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Again three main regions occur (Fig. 3); in region I the perturbation 
potential is zero and regions II and III correspond to the orifice being 
intersected by, and lying wholly within, the forecone of P respectively; 
further (mathematical) subdivision is occasioned in evaluating the 
potential in region II, and the ellipse x?/B?+-r? = a* and the hyperbola 
x?/82—y? — a® bound these subdivisions; as would be expected, the final 


Fic. 3. Integration regions and sub-regions for an axially 
symmetrical nozzle flow. 


results obtained are the same for the whole of region II. The actual 
process of evaluating the potential is very similar in each of the 
regions, and a brief outline of the steps involved will be given for 
region ITI. 

Taking the potential as the difference between the potentials due to 
the complete plane and to the orifice, we have 


r+z/B xX 


(n?+¢?—a?) dy 


a A 
df 
(5.3 


where X = A = 


The first of the repeated integrations may be reduced to standard forms 


I 
é 
é 
i 


| in 
I 
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in each case and the first term integrated completely, whence 


0 

— (7 + 
+r 
(5.4) 

It is convenient to obtain the velocity components at this stage since 
the potential itself is not required: thus 


+ 
a*—»?)+r 

+{rsin-4 +rsin- dyn. (5.5) 

It may be mentioned in passing that the limits of integration in 

Region II involve both x and r, but, in differentiating under the 

integral sign, the terms obtained from differentiating the limits dis- 

appear. In the evaluation of the integral occurring in (5.5), the terms 

are considered in pairs as bracketed, and the substitution 7 = 2at,/(1—t?*) 

is used. The range 0 < » <a corresponds to either 0 <¢ < 1/¥2 or 

1 >t > 1/v2, and, by using a different range for each term in a pair, 

the pairs each reduce to a single term with limits 0 and 1, and it is not 
difficult to show that 


B 
ra2 2 
32 (om) — Blom} — — a+r) Km) + 
(A+r 2a 2a 
where g= m= = — 


and K, HE, and II are complete elliptic integrals of the first, second, and 
third kinds respectively (as defined by Legendre); in Region III, 
0 <m <1 and A is imaginary. The expression given for @¢/ér may 
be reduced considerably, but it is quoted in the form in which it is 


| 
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obtained from the steps listed. Since complete elliptic integrals of the 
third kind may be expressed in terms of integrals of the first and second 
kinds, (5.6) can be materially simplified. 

The reduction is most conveniently performed using the Jacobian 
addition theorem in the Legendre notation; effectively this is a 
relation between three elliptic integrals of the third kind having similar . 
moduli and different parameters, and it is found that, if the integrals 
in (5.6) constitute two terms (with complex parameters), the coefficients 
occurring in the theorem are just those required for the reduction, and 
that the third integral has a real parameter. Since this integral is also 
complete, it may be expressed in terms of integrals of the first two kinds 
using the relations given by Legendre in his T'raité des Fonctions 
Elliptiques. 

By working in terms of elliptic functions, i.e. by expressing a 
parameter n as —k?sn*a (k being the modulus), the question as to | 
whether the parameter is real or complex is avoided since Richelot has 
shown that any complex number may be expressed in such a form (with 
a complex). The evaluation of the perturbation potentials is rendered 
quite straightforward by following this procedure although the algebra 
is rather lengthy, and the results are quoted: 


region I: ~ = 0; 
region II: 
+41 4a) BW, a) 
region ITT: 
= 


4K(m){2 2r? [2 


E(m', 


{K(m)—Bom) + (Fe rom’), 


2x 


Ng B? 
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where 2 
m= ,/(1 ma) = @ = vy = sin 


and F is an incomplete elliptic integral of the first kind. 


4293-4604 (9042) 


506— 745— 1018 1320—-165!—2000-2364- 2720-3083 365-3945 (6740) 


“O— 24— 449— — — (7942) 


(i009 (72 (1724) (2511) 385) = (6589) 


Fic. 4. Radial velocity component in a specific axially symmetrical nozzle 
flow at a Mach number M = v2. (Values of é¢/éx in brackets.) 


Fic. 5. Divergent and convergent axially symmetrical nozzle 
flows corresponding to the results of Fig. 4. 


The value of @¢/@x at r = 0 reduces to an expression which can be 
obtained otherwise; the values of é¢/ér for regions II and III agree on 
their common boundary, and the radial perturbation velocity is zero 
along the axis and on the boundary between regions I and II, i.e. the 
transition between the initial stream and the modified flow is smooth. 
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The results have been computed for a Mach number M = v2 using 
the Smithsonian Tables, and are given to a convenient scale in Fig. 4. 
The two flow patterns corresponding to this distribution are shown in 
Fig. 5; a multiplying factor 10-*U has been used in conjunction with 
the results of Fig. 4, and thus it can be seen that the assumptions of the 
linearized theory are not satisfied over the whole flow patterns given. 


6. Conclusion 


The choice of a density function is limited by the proviso that the 
associated perturbation velocities must conform to the restrictions im- 
posed by the linearization. It is found that, for a uniformly dense source 
plane containing a circular orifice, infinities exist in the perturbation 
velocities, and consideration of a generalized distribution suggests that 
a sufficient condition to avoid similar infinities is that the source density 
should be zero on the edge of the orifice; this condition holds for the 
example given. 

Investigations made in connexion with aerofoil theory suggest that 
a smooth transition between the undisturbed and disturbed regions can 
be expected for density functions having a finite gradient at the orifice 
edge, and that an infinite gradient, e.g. f(y, ¢) = corre- 
sponds to a discontinuity in the velocity vector at this boundary 
although it may well transpire that further restrictions must be 
imposed on the choice of a suitable function. It is remarked that the 
ultimate nature of the flow pattern is similar to that produced by a 
plane of similar density but without the orifice. 

The method is not confined to axially symmetrical flows although 
naturally the algebra will become considerably more involved for flows 
not possessing axial symmetry. 


7. Virtual source distributions in aerofoil theory 


The solution of the linearized equation for three-dimensional aerofoil 
problems has been considered by several writers, and the application 
of Green’s theorem to such problems (3) involves the use of equivalent- 
layer distributions of sources and doublets; the results given in this 
section are derived from true virtual distributions. 

If the initial stream of (supersonic) velocity U is incident parallel to 
the z-axis with y spanwise and z vertical, it is within the limits of the 
approximations made to assume that a slender wing at small incidence 
occupies a plane z = constant as far as evaluating the perturbation 
velocities is concerned. In using equivalent layer distributions, the 


ON SUPERSONIC FLOWS 301 


distribution is taken to replace the area occupied by the wing. If this 
distribution is taken to be virtual, any plane can be taken to contain 
the wing, and a series of wings can be associated with each distribution. 

The wing examples considered below have leading edges of the super- 
sonic type, i.e. the leading edge lies at all points ahead of the local Mach 
cone, and hence any discontinuities in the perturbation velocities on 
crossing the wing can be supported by the wing; as in the equivalent- 
layer theory the flows above and below such wings are considered 
separately. 

A simple example derived from a line source is given to clarify the 
method. 


8. A specific virtual line source distribution 

Consider a line source of density f(n,) = ./(¢?—a?) disposed along 
the negative axis of { (or z) with —o <{< —a. The potential at 
P (x,y, 2) is 


V(—a*) de 


—y?) 
the lower limit being obtained as the (negative) intersection of the 
forecone of P with the C-axis. 

The associated wing lying in the plane z = 0 will have as its leading 
edge the curve in which the plane cuts the cone envelope of the line 
source, i.e. the hyperbola x?/8?—y? = a*; the potential (8.1) is assumed 
to hold for the region lying above the wing. 


Cdl 


where s? = 2?/8?—y?, whence 


This line source, therefore, is associated with a flat inclined plate with 
a hyperbolic leading edge: 


and, taking k = a/s (0 <k < 1), we have 


= 


3 
4 7 
= 
= 
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where L’ is a complete elliptic integral of the second kind with the 
complementary modulus k’. 

The perturbation velocities occurring on the upper surface of a flat 
plate at incidence a and with leading edge 


> 
(x > fa) 
would be 


The flow pattern below this plate can be inferred immediately from 
these results, but it could be obtained from a similar line source of 
opposite sign disposed along the positive axis for +a <{[<+o, 
although the flow would not occur in the simultaneous presence of 
both distributions. 


9. A symmetrical delta wing at incidence 

The problem of a given symmetrical wing with a supersonic leading 
edge may be approached directly from the equivalent-layer standpoint, 
but the inverse mode of approach, by means of which the nature of a 
wing corresponding to a postulated distribution is determined, is useful 
since exact solutions for winglike surfaces may be obtained using 
comparatively simple mathematics. 

If a ‘triangular’ virtual source distribution is located in the lower 
half (7, £)-plane with vertex at the origin and sides equally inclined at 
an angle 7 to the ¢-axis (Fig. 6), the leading edge of the associated wing 
in the plane z = 0 is of the delta type given by x = +fycosr, x > 0. 
The potential is evaluated in two regions depending on whether the 
forecone of a point intersects the source triangle on one or on both sides, 
and these two cases correspond to the regions lying ahead of and within 
the Mach cone of the vertex respectively. It is found that a uniformly 
dense distribution corresponds to a flat inclined delta wing, and it will 
be recalled that Puckett (4) has considered this case using a uniformly 
dense equivalent-layer distribution. The trailing edge may be selected 
quite arbitrarily without affecting the validity of the solution provided 
it is of the supersonic type. 

In considering a wing of finite thickness, e.g. a symmetrical wing, this 
is no longer true, and the leading and trailing edges and the camber 
cannot be selected independently. A symmetrical wing at zero incidence 
is a slender surface with sharp edges lying in the plane z = 0 (say). 
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A solution of the potential equation determines the downwash on the 
wing, and the wing surface is given by 

w w 1 

— = . 9.1 

The solution of this equation admits of the addition of an arbitrary 
function of y which is determined by the condition that the wing 


A 


Fig. 6. Integration region for a triangular plane 
source distribution normal to the stream. 


surface intersects the plane z = 0 along the leading edge. Since the 
camber and leading edge are already specified at this stage, the choice 
of trailing edge is restricted to the curve or curves in which the surface 
corresponding to (9.1) intersects the plane z = 0 again. 

If a source distribution function 


f(y, 6) = AP + Br?+C 

is taken for the triangle cited, it is found that the upper surface of a 
symmetrical wing at zero incidence associated with this distribution 
for the plane z = 0 is 

Uz = —(x+By cos 7){aa(a+ By cos 7)+bB?y?—ac?}, (9.2) 
where c is the centre-line chord and a and 6 are linear functions of A 
and B; the upper sign refers to positive values of y. Since b is a function 
of the arbitrary coefficients A and B, it may be expressed as 


b = a cos?7, (9.3) 


P 
— 

U 
y 
SS 
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from which it is seen that, since 8 and 7 always occur in the combination 
Bcosr, a value of 7 can always be found for a given f to preserve the 
geometrical shape of the wing. 
The semi-angle w of the delta wing is related to 8 and +r by 
cot w = Bcosr. 

The first bracket in (9.2) equated to zero gives the leading edge, and 
the second the possible trailing edges as a one-parameter family of 
curves. The nature of these trailing edges is 
shown in Fig. 7; to preserve the validity of 
the velocity components over the whole wing 
surface, the trailing edges must be wholly 


F supersonic, and the curves shown are all of 
this type. 
RECTILINEAR The lift and drag coefficients for a sym- 


TRAILING EDGE etrical wing at incidence are 


4 
Cop (w, U;) ds, 
S 


where suffixes i and ¢ refer to the perturba- 
hs, tion velocities occurring on the upper surfaces 
, of a flat plate with the same plan-form and 
incidence as the wing, and on the non- 
incident wing respectively. 

The lift coefficient is therefore dependent 
Fic. 7. The family of super- 5) incidence alone, and the drag coefficient 
sonictrailingedgesforaspeci- 
fic symmetrical wing with iS the sum of the coefficients due to the flat 

delta-type leading edges. plate at incidence and to the profile at zero 

incidence. 

Taking A = }, the expressions for the lift and drag coefficients may 
be integrated completely; the leading and trailing edges of this par- 
ticular (non-incident) wing are, for y > 0, 

x = Bycosr, x+4By cost = c, 
and the perturbation velocities occurring on the upper surface of the 


wing are w = a(c?-+ cos*r— 322), 


{x?(6—53 cos*r + 33 cost) —aBy($-+ 12 


B?y?(83 cos*z7 — $)cos*7 — 2c? sin‘7}, 


= 
sin5r 


on 
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= cos?r)(cos 7 sin tT) — 


+tan- {a?(6—3 cos*r-+- 3} cos*r) + 

cos*r — 3)cos*x— 2c? sin‘r}], 
where suffixes a and } refer to the regions lying ahead of and within 
the Mach cone of the vertex respectively. 


395 — 
390 VA 330 
385 3:25 
3:80 7 320 


375 3-15 
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1300 


56077020 30 40 7 00 60.10 80 90 
Fic. 8. Incidence and profile drag-functions for a specific symmetrical wing. 


The lift and drag cov:.icients are found to be 
2sinr 
Brsinz\./(1—} cos?) 


62 a2 


a 
D? 


cos~1(4 cos7)—7] = g(r) 


Cr 


Cp = f(r) 2 +9(7) 


where « is the incidence, @ is the thickness ratio defined as the maximum 

thickness in the centre-line section divided by the chord, and f(r) and 

g(r) are functions of 7 given graphically in Fig. 8; the algebra leading 
3695.2.2 x 
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to the determination of f(r) is lengthy and the function is complicated, 
involving factors containing polynomials of the twelfth degree in cos r, 
but it is a check on the accuracy that the value at 7 = 0 is finite since 
it is evaluated as the limiting value of a function h(r)cosec*r as + > 0. 

From Fig. 8 it is seen that the least profile drag for a given thickness 
ratio is experienced by a wing having the greatest semi-vertical angle 


whereas the incidence drag is least for a wing whose leading edge lies 
in the Mach cone of the vertex. 


SECTION 
W=0 


CENTRELINE SECTION 


Fic. 9. A symmetrical delta wing at zero incidence showing 
the centre-line and maximum thickness sections. 


The geometry of the wing with semi-vertical angle w = cot-!0-5 at 
a Mach number M = 2 (7 = 47) is given in Fig. 9. The section w = 0 
is centrally symmetrical, and from the quadrant given it is seen that 
the wing has a sharp ridge along the longitudinal axis. 

A special case of the above is the flat plate at incidence (9 = 0), and 
the results obtained by Puckett (4) for the velocity components are in 
agreement with the values obtained from the virtual distribution. 


10. Remarks 


The use of virtual distributions has enabled the source-sink method 
of supersonic flows to be extended to the consideration of nozzle-flow 
problems. The results obtained for the specific axially symmetrical 
flow represent the simplest available exact solution of the linearized 
equation for this type of problem, approached from the virtual distribu- 


tion standpoint; the method is quite general and not confined to flows 
possessing axial symmetry. 
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Certain types of aerofoil problem are amenable to solution by the 
method, which includes the equivalent-layer distributions as special 


cases. 
I should like to thank Professors H. Levy and W. G. Bickley for their 


interest and encouragement; my very sincere thanks are due to Mr. C. H. 
Kebby for many interesting discussions and for his generous help. 
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NON-HOMOGENEOUS LINEAR FORMS 
ASSOCIATED WITH ALGEBRAIC FIELDS 


By L. E. CLARKE (London) 
[Received 15 November 1950] 
1. Let L,,..., L,, be any n linear forms in n variables 2,,..., x,,, with real 
coefficients and determinant A, not zero. It is known that there exists 


a constant A,,, depending only on n, such that for any n real numbers 
the inequality 


| < A, |A| (1) 
is soluble in integers 2,,..., z,. A well-known conjecture of Minkowski 
(proved as yet only for n < 4) asserts that this is true when A,, = 2-". 


Tschebotareff has provedy that it is true for any constant A, which is 
greater than (v2)-”. 

Now consider the more general situation where the forms J,...., L,, 
consist of r forms with real coefficients, and s pairs of forms with con- 
jugate complex coefficients, where r+-2s = n. We now allow 4y...., a, 

to be real or complex numbers, stipulating that those corresponding to 
real forms shall be real, and those corresponding to conjugate complex 
forms shall be conjugate complex. As one sees already from the case 
r = 0,8 = 1, no result exactly analogous to (1) can be valid. We have 
only to take the forms 

L, = L, = «,—i0x, (0 > 0), (2) 
and the numbers a, = 3(1+76), a, = 3(1—7@), to get an immediate 
contradiction to (1). For then the product on the left is always at least 
3(1+6?), and |A| = 26, giving a contradiction if @ is sufficiently large 
or sufficiently small. 

The simplest sets of linear forms which suggest themselves in this 
connexion are those associated with algebraic number-fields. Let # be an 
algebraic number of degree n, with conjugates 0™,..., Let 0%,..., 
be real, and the others complex, and let the latter be numbered so that 
and &°+*) are complex conjugates, for v = r+1.,..., r-+s (r+2s = n). 
Let wj,..., w, be a basis for the algebraic integers of the field k(@), and 
let w,..., w, for 7 = 1,...,, be the corresponding basis of the field 
k(@). The linear forms 

L; = (j = 1,..., ) (3) 

+ The proof is given in Hardy and Wright, Introduction to the Theory of Numbers 
2nd ed. (Oxford, 1945) 396-8. 

Quart. J. Math. Oxford (2), 2 (1951), 308-15. 


ON NON-HOMOGENEOUS LINEAR FORMS 309 


consist of r real forms and s pairs of conjugate complex forms. As 
2y,..., &, take all rational integral values, L, represents the general 
integer of k(O™) and Ly,..., L, are the algebraic conjugates of L,. The 
determinant of the linear forms is + vd, where d is the discriminant of 
the field k(@). 

The object of this paper is to establish some results which are partially 
analogous to (1), when J,,..., L,, are linear forms of the kind just defined, 
and 4,,..., @, are as stated earlier. Such results, as we have seen, cannot 
hold with |A| or ,/\d| on the right, but I shall prove that they do hold 
when this is replaced by a suitable power of |d|. The simplest result 
of this kind is Theorem 1, which asserts that the inequality 


< || (4) 


is soluble in rational integers 2,,..., x,, where , is a suitable constant 
depending only on n. 

The example (2) considered above, with 6 = vm, where m is a large 
positive integer, square-free and congruent to 1 or 2 (mod 4), shows 
that the exponent $n in (4) cannot be improved upon when n = 2, 
r = 0,8 = 1. It is very probable that a better exponent will suffice in 
all other cases. An improvement is obtained in Theorem 2 for all cases 
in which n is a prime greater than 2, and here the exponent $n is 


replaced by jn—(n—2)/(n—1). 


Recent work of Davenport} has thrown light on this question in the 
three cases when n = 2,r = 2,8 = 0,orn = 3,r = 1,8 = l, orn = 4, 
r = 0, s = 2. In this work it has been proved that an inequality 
analogous to (4), with the exponent $, cannot be valid in any field of 
these three types, if the constant p,, is sufficiently small. 

The simplest type of field for which the best possible exponent re- 
mains unknown is that for which n = 3, r = 1, s = 1, i.e. cubic fields 
of negative discriminant. The simplest such fields are those generated 
by m*, where m is an integer, not a perfect cube. To avoid complications, 
I shall suppose that m is square-free, and that 1, 0, 6? is an integral 
basis for the field k(@), where 6 = m+. Under these conditions, I prove 
in Theorem 3 that the inequality is soluble with a constant multiple of 
\d\§ on the right-hand side. Further, in Theorem 4 I prove that there 
are infinitely many fields k(m*) for which any result with a smaller 
exponent would not be true. 


+ Acta Math. 84 (1950), 159-79, and Trans. American Math. Soc. 68 (1950), 
508-32. 
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2. Let L,,..., L, be the linear forms associated with an algebraic 
number-field of degree n and type (r, 8), as defined in § 1, so that 
L,,..., L, are real and L; is the complex conjugate form to L,,, for 
j= r+s. 

Ciena 1. There exists a number p,, depending only on n, with the 
following property: if ay,..., 4, are any numbers with ay,..., a, real and a, 
the complex conjugate of a;,, for j = r+-1,..., r+-s, then the inequality (4) 
is soluble in rational integers x,,..., X,. 

Proof. By the restrictions imposed on 4j,..., a,, there exist real 
numbers 2},..., 2;, for which 

... L, = —a,. 
On changing ,,..., 2, into x,+2},..., %,+2},, we see that it suffices to 
prove that the inequality 


is soluble in real numbers 2,..., x, satisfying 
=2,(modl), ..., 2x, =2;,(mod1). (6) 


Following Hermite, we associate with the linear forms J,,..., L,, the 

positive definite quadratic form 
Q(x,,...,%,) = L3+...4 124+ |L,4,|?+...+ |*. (7) 
To calculate the determinant of this form, we write for the moment 
= = (¥;+1Z;) (j | 8), 
where Y;, Z; are real linear forms. The determinant of the n real forms 
Y,... ., Y,, Z,,..., Z, is obtained from the determinant of Z,,..., L,, 
by multiplying by 7°, sal so is +i*vd. Since 

the determinant of Q (being obviously positive) is simply |d|. We note 

We apply an integral unimodular transformation to the variables 
X15... L, 80 as to reduce the form Q(2,,...,2,,) in the sense of Minkowski. 
After such a transformation, we can write 


n n 
p=1 q=1 
wheret Ogg S K (9) 
and < K,|d|, (10) 


for some constant K,, depending only on n. The effect of an integral 


+ Minkowski, Ges. Abh. ii. 53-100. 
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unimodular transformation on the conditions (6) is to transform them 
into conditions of the same kind on the new variables. 
We now express Q as a sum of squares, in the form 


+5, +... +8, 23. 


By comparison of coefficients, we have 


By < Bn < Xnn- 
We choose 2,,, successively, to satisfy the conditions (6) and 
also to satisfy <4, 


+4 < 3, 


We then ensure that (6) holds, and 
Q(X HBy+--- < (11) 


In order to obtain an upper bound for «,,,, in terms of |d|, we must 
use the fact that the linear forms J,,..., L,, correspond to an algebraic 
number-field of degree n. For all rational integral values of 2,,..., %», 
the product L,...L,, is the norm of the algebraic integer L,, and so is 
a rational integer, and is not zero unless 2,..., 2, are all zero. Hence, 


by (8), we have Q(2x,,..,2,) 


for all integers 2,,..., X,, not all zero. In particular, «,, >”. Hence, 
by (9) and (10), Onn < K,,|d|n3-. (12) 


Returning to (11), we have found real numbeis 2,..., 2, which satisfy 
By (8), this implies 
|Ly...L,| < |d| 
which proves (5). 
3. THEOREM 2. For each prime value of n, there exists a number ,, such 
that, with the hypotheses of Theorem 1, the inequality 


< |d|” (13) 
is soluble in integers 2,,..., ,, where 
v = 4n—(n—2)/(n—1). (14) 


Proof. The improved result depends upon obtaining, for prime values 
of n, a better estimate for a,,, than that given in (12). The rest of the 
proof is the same as before. 


r 
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The minimum of Q(2,,...,x,) for rational integers 2,,..., %,, i8 a,,, and 
this minimum is obviously n, and is attained for those integers f,..., «* 
which make L, = 1. By Minkowski’s definition of a reduced form, a,, 
is the minimum of Q for rational integers 2,,..., 2, which are not propor- 
tional to af,..., For such values of x,,..., the number is an 
algebraic integer of k(@™), but is not rational. Since n is a prime, L, 
must be an algebraic number of degree n, and so 


I] (L,—L;) # ©. 


h<j 
It is easily seen that this product is a rational integral multiple of vd. 
In fact, by (3), we can express the powers of L; in the form 
= Way +... 
forj = 1,...,nandh = l.,...,, where the 2’s are rational integers. Hence 
= || x |larjall, 
whence [] |L,—L;| = xv|d|, where x is a rational integer. 


h<j 
For any such values of 2,,..., X,, we have 


Q(a,,...,%,) = 124+ |L,,,|?+...+- |? 
1 


h<j 


1 
2, (|Z;,|+- |Z;|)? 


> TT 
h<j 
In particular, > 


and so (9) and (10) give 
n 


Onn 


Using this instead of (12), we obtain (13) with the exponent v given in (14). 
4. Let m be a positive square-free integer greater than 1, and let 
6 = m*. We assume that 1, 0, 6? form a basis for the integers of the 
field k(@), as is the case} if m = 0 (mod 3) or m = +2 or +4 (mod 9). 
The discriminant of the field is —27m?. Let L, be the linear form 
L, = 6y+ (15) 
and let L,, L, be the algebraically conjugate forms. 
+ See Sommer, Vorleswngen tiber Zahlentheorie (Berlin, 1907) 261-2. 
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THEOREM 3. If a, is any real number, and ay, a, any pair of conjugate 

complex numbers, then there exist rational integers x, y, z such that 
< mi. (16) 
Proof. As before, it suffices to find real numbers 2, y, z which satisfy 
conditions of the form 
x = 2’ (mod 1), y =y’ (mod 1), z = 2’ (modl), (17) 
and which satisfy also |L, L, L,| < mi. (18) 
We write 
X = «4+ by+- Y = Z =z. 
Then L, = X, and 
L, L, = 6y?+ 042?— — Oxy 
= X*—30XY 
We choose z, y, x in succession to satisfy (17) and 
IZ}<4 |X| <}. 
Then Ly Ly| < 4(1-+30-+ 362+ 268), 
Since 6? > 2, this is less than 64 = mf, and the result (18) follows. 

It may be observed that the hypothesis that m is square-free is not 
essential, provided that m is cube-free and satisfies m = 0 (mod 3) or 
m = +2 or +4 (mod9). In this case, if we write m = m,m}, where m, 
is square-free, a basis for the field k(m*) is givent by 1, 0, 63, 63 6,, where 
6, = mj and 6, = m}. The discriminant of the field is now given by 
d = —27m? m3. A similar calculation shows that the conclusion (16) is 
still valid, apart from a constant factor on the right-hand side. This 
applies equally whether we take mi or di on the right. 

5. We now consider the same linear forms as in § 4, but assume that 
m is of the special form /*+-1, where / is a large even positive integer. 
We postpone until § 6 the consideration of whether such values of m have 
the properties postulated at the beginning of § 4. We simply consider 
the linear form defined by (15), and its conjugates, and prove that the 
exponent # on the right of (16) is the best possibie. 

THEOREM 4. Let m = [3+-1, where l is a sufficiently large positive even 
integer, let @ = m+, and let L, be the linear form defined by (15), and L, 
and Lz be its algebraic conjugates. Suppose that C < 3;. Then there exist 
a real number a, and a pair of conjugate complex conjugate numbers az, C3 


such that > (19) 


for all rational integers x, y, 2. 
t+ See Sommer, loc. cit. 
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Proof. We choose a, = 4(1+6-+6?), etc., so that the conditions (17) 
become x = (mod 1), y = (mod 1), z = 4(mod1). We have there- 


fore to prove that \L, Ly > 8Cm* (20) 


for all odd rational integers x, y, z. 

We note that (@—1) = 72+10+ 6, 
and consequently @—/ is a unit, since its reciprocal is an integer. 

We define an integral unimodular transformation from variables z, 
y, z to variables x’, y’, z’ by 

a’ + By’ + 62’ = (8—I)(x+ by+ 62), 
which is equivalent to 
= —lx+(+1)z, y = «—ly, 2’ = y—lIz. 
Since / is even, we have 
x’ = z(mod 2), y’ = x (mod 2), 2’ = y (mod 2), 
and so x’, y’, 2’ are all odd if and only if x, y, z are all odd. This trans- 
formation changes L, into L,, where 
L, = 

The transformation does not change the value of L, L, L, since @—1 is 


a unit of the field k(@). 
Suppose that there are odd integers x, y, z which violate (20), so that 


|L, L, < P, P = 8Cmt‘. (21) 
After applying the above transformation, or its inverse, a suitable num- 
ber of times, we obtain odd integers x, y, z which satisfy both (21) and 
< |L,| < }. 

Then |\L, L,| < 3P/(@—1l) = 3P(l?+10+-6?) < 9P6?. 
Thus |L,| = < 30VP. (22) 

We now express y and z in terms of L,, L,, L, by the equations 

L, = «+ L, = L, = wOy+w6z. 


30% = 
We now write L, in the form 
L, = 
Since |L,| < 4 and x+ly+/*z is an odd integer, we have 
> §, 
i.e. ly+(0+l)z| > = 6?) > (24) 


7) 
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Since P = O(/*) and 6 = O(l), where the O-notation refers to / > 00, 

it follows from (22) that |Z,| = O(/°), and then from (23) that 
y=O0), z=O0(. 
Hence (24) can be written as 
21? < = \y+26z|+O(I-). 
Now, by (23), 
30(y+20z) = 3L,+(w?+ 2w)L,+ (w+ 2w?) Ls. 
Since |w?+2w| = |w+2w?| = v3, 
we obtain from (22) 
30|\y+202| < 14-66,/(3P). 

Hence 2? < 2,/(3P)+O(l-). 
Since P = 8Cm* = 80C(/°+-1), this is impossible for large 1 if C < 3. 

6. In order to prove that the result of Theorem 4 is relevant to that 
of Theorem 3, and thereby to prove that the exponent $ in Theorem 3 
is the best possible, we need to know that there are infinitely many even 
integers / for which /*+-1 is square-free, and = 0 (mod 3), or +2 (mod 9), 
or +4(mod9). The existence of such numbers is easily established as 
follows. We take / = 6t+-4; then / is even, and 1 = 2 (mod 9). Also 
B4+-1 is not divisible by 2 or 3. The factors /+-1 and /?—1+-1 of B+1 
are relatively prime, and it suffices to find values of ¢ for which each 
of them is square-free. If t < 7’, any prime p > 5 for which p? divides 
1+1 or /?—Il+1 must satisfy p = O(7). Also 1, and therefore also ¢, 
must be congruent to one of at most three numbers to the modulus p?. 
Hence the number of values of t < 7' for which either /+-1 or /?—1-+-1 
is not square-free is less than 

summed over primes p > 5 with p = O(7'). Since the number of primes 
p in the sum is O(7'/log 7’), and since 
> 3/p? <1, 
pos 


the existence of infinitely many / with the desired property follows. 
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ELECTRICAL NOTES 


By F. B. PIDDUCK (Keswick) 


XIII. DIFFRACTION OF LIGHT BY A SEMI-TRANSPARENT 
SHEET 


[Received 27 December 1950] 


THE elementary theory of diffraction, based on Fresnel’s zones, was 
developed by Kirchhoff so as to give the intensity of light at any point 
due to an aperture of any form in an opaque sheet, in terms of an integral 
taken over the aperture. Though the perfected zone theory is not rigor- 
ous, it solves more problems in its own field than the electromagnetic 
theory, which depends on the solution of partial differential equations 
with boundary conditions. That field is discussed below. Kirchhoff’s 
formula was derived as a first approximation from the electromagnetic 
theory by L. de Broglie (1). Magnus (2) showed that Sommerfeld’s 
problem of a straight edge can be solved by first calculating the currents 
in the sheet. If the linear dimensions of the sheet are large compared 
with the wave-length, the distribution of dipoles can be taken to be that 
in an infinite sheet of the same material and thickness. Kirchhoff’s 
formula can be established in this way, and extended to screens of any 
degree of transparency (3). The screen may in particular be a glass 
plate (4). 


General characteristics 

The fundamental problem is that of a thin disk of homogeneous con- 
ducting dielectric whose linear dimensions are small, but much larger 
than the wave-length A. The source P, emits a wave in which the electric 
intensity at distance r is A exp(twt—vkr)/r, where k = w/c = 27/A. Let 
P’ be a point on the flat surface of the disk nearest to Pj, and P a point 
of observation such that either P,, P’, and P are nearly in a straight line 
(transmission pattern) or such that P, P’, PP’ make nearly equal angles 
with the normal at P’ (reflexion pattern). We consider in the first 
instance light polarized in and perpendicular to the plane of incidence, 
from which other polarizations can be obtained by superposition. The 
electric intensities at P due to the moments induced in the disk are then 
found, by evaluating expressions given in (3), to be 


E’p> = —(1—T)KA _ in the transmission pattern, (1) 
E> = RKA in the reflection pattern, ( 


bo 


Quart. J. Math. Oxford (2), 2 (1951), 316-20. 
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where exp(—ikv,—ikr) dS. (3) 
0 


Here ¢ is the angle of incidence, » = P, P’, r= P’P, dS an element 
of area of the disk at P’, and integration is over the disk. 7’ and R are 
the fractions of electric intensity transmitted and reflected by an infinite 
sheet of the same material and thickness. It is here that the specific 
properties of the material make their appearance. For a perfectly con- 
ducting sheet 7’ = 0 and R = 1; but in general 7 and R are different 
for different polarizations, and the patterns are of different absolute 
intensity, though the maxima and minima are in the same position. 
The transmission pattern is the same for all thin opaque sheets, of what- 
ever material they may be composed. 

The factors 1—7' and R have been taken out of the integral in equa- 
tions (1) and (2). We shall find how far this is legitimate for a disk 2 mm. 
in diameter and of thickness b = 0-1 mm. when the electric intensity is 
in the plane of incidence. If n is the (complex) refractive index and ¢’ 
the (complex) angle of refraction, 

Tl — 4n cos ¢ cos ¢’ exp(ikb cos ¢) 

(n cos cos exp(inkb cos ¢’)— 
—(n cos 6—cos exp(—inkb cos ¢’) 
cos ¢’)—exp(—inkb cos ¢’) | 
(ncos¢-+cos ¢’)? exp(inkb cos 
—(ncos d—cos $’)? exp(—inkb cos ¢’) 


> 


We write n = Ny—W, (4) 
where m, and y are real and positive. The cases of no appreciable error 
are: 

Fraunhofer diffraction. T and R are constant and (1) and (2) are 
correct. 

Metal sheets. v lies between 2-5 and 5 for the common metals, and 
kb = 1070 for sodium light. 7’ is extremely small and 


R = (ncos¢—cos ¢’)/(n cos cos 9’), 


which varies only slowly with ¢. 

Normal incidence. Consider in particular a transparent sheet of refrac- 
tive index 1-5. If ¢ is the greatest angle of incidence, and the source is 
x mm. from the sheet, kbcos¢ varies from centre to edge by approxi- 
mately kbd?/2 = kb/2x*. The variation of nkbcos¢’ is 1/nth of this. If 
the greatest of these variations is (say) 0-05, 7'is approximately constant. 
This makes x = 103. If the sheet is exactly semi-transparent, exp(vkb) 
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is of the order of V2, and exp(vkb cos ¢’) varies very little. Equation (1) 
is accurate if the source is at least 10 cm. away. 

The sheet has been taken in this theory to be uniform in thickness and 
refractive index, and the light homogeneous. A small change of ¢ then 
makes the whole pattern change in intensity, without much effect on the 
positions of maxima or minima. The sheet, that is, is an interferometer, 
and the diffraction pattern follows suit. Accurate figuring is required to 
reveal this, and in practice no rapid change of intensity is noticed. 


Straight edge 
Let a semi-transparent sheet occupy the part x < 0 ofthe planez = 0, 
and let a plane wave in which the electric intensity is exp(twt+-tkz) fall 
on it normally from above. Let P be the point (x7, 0, —a), at distance x 
from the geometrical shadow and a below the sheet, where x is much 
smaller than a. We can, as usual, take Kirchhoff’s integral over the 
whole half-plane, with r inside the integral. A/r, is replaced by 1. The 
electric intensity at P due to the sheet is thus 
—i(1—T )expiwt exp(—ikr) dS 
A | r 
The calculation can be put very concisely by using the properties of 
Bessel functions. Let P’ be the point (x—z’,y’,0). Then dS = dz'dy’ 
andr = (R*+y’'*)!, where R = (a?+<2’2)t. The result of integrating with 
respect to y’ is —iwH?(kR). If kR is large, 
HP(kR) Rt, 
so that Re 
—(1—T)(1+ exp iwt | exp(—ikR) dx’ 
Aki Ri 


= 


Most of the integral comes from the smaller values of x’, so that we can 
replace R by a+-2’*/2a in the numerator and by a in the denominator. 
Changing the variable from x’ to k+x’/2a*t, we have 

E, = 
where 


Me 2 S== | singtds. 6) 


Adding in the incident wave exp(iwt—ika) we have for the total electric 
intensity at P 


E, = 1+ 7+(1—T7)(1+7)(C—iS)]. (6) 


), 


= 
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The transmission factor 7’ is complex in general, the modulus deter- 
mining the intensity of the transmitted light and the argument its phase. 
Whatever it is, the intensity of light in the diffraction pattern can be 
found by a suitable choice of the base point in Cornu’s spiral. In that 
curve S is plotted against C, and the values of € are marked on the curve. 
The curve has two winding points A and B, corresponding to = —oo 
and oo. If P is any point on the curve, the intensity of the light dif- 
fracted by an opaque screen is proportional to A P?. For a screen of any 
transmission factor 7’, the intensity is, from equation (6), 


(s—S)*], 
where c and s are real numbers such that 
c—is = —(14- T)(144)(1—7). (7) 
Thus, if Gis a point whose coordinates are c and s in the plane of Cornu’s 
spiral, the intensity of light is proportional to GP?. The pattern has 
fringes on both sides of the geometrical shadow, except in the cases 
T = 0 (opaque screen), where G coincides with A, 
T = 1 (no screen at all), where G is at infinity. 


The diffraction pattern of a semi-transparent sheet with a straight 
edge, photographed by Zernike (5), is accounted for by this theory, in 
so far as the material of the sheet may be considered as a homogeneous 
conducting dielectric. I have obtained good patterns by reflection from 
the face of a prism, using an ‘Osira’ sodium lamp and a panchromatic 
film, with an exposure of 10 minutes. They are hardly distinguishable 
from the transmission pattern of an opaque screen with a straight edge. 


Forms of Babinet’s principle 

Let a small disk D be cut out of an infinite sheet, and let S be the rest 
of the sheet. When D and S are used separately as diffracting screens, 
let E be the electric intensity of the incident light and E’p, Eg those of 
the moments induced by it in D and S. The total electric intensities in 
the diffraction patterns are Ep = E+E’ and Eg = E+Eg. By our 
approximate hypothesis, the two sets of moments together make up the 
moment in the uncut sheet. Thus E, = E+E‘+Eg, where Ez is the 
electric intensity transmitted by the uncut sheet. From these equations 


= E+E, (8) 
and E, = E,—E>. (9) 
In Fraunhofer’s diffraction phenomena E = E, = 0 except at the cen- 
tral spot, and Ep, = —Eg. The intensity of the diffraction patterns of 
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the disk and the rest of the sheet is the same except at the central spot, 
whatever the transparency. Equation (9) gives the electric intensity in 
the pattern of a perforated sheet when the effect of the complementary 
disk has been calculated from (1) and (3). 

For reflected patterns, with the point of observation screened from 
the direct light, Ep = E’p, Eg = Es, and 

Ep+Eg = Ep. (10) 

Scope of the theory 

Equations (1) and (2) hold only when the angle of diffraction is small. 
The density of electric moment is constant to within a few wave-lengths 
of the edge, and this constant part contributes most of the electric inten- 
sity of the sheet at small angles. The edge, on the other hand, contri- 
butes a finite proportion at large angles. This is why Kirchhoff’s formula 
fails at large angles. It is, perhaps, less well known that it fails for small 
apertures. This was first pointed out by Rayleigh (6), (7), who cited 
an observation of Fizeau that a very fine slit polarizes the light that 
it transmits. The effect is easily seen with a sodium lamp and a pola- 
roid disk. For these reasons Kirchhoff’s formula cannot be used with 
gratings, in which the angle is large and the width of the slit often less 
than a wave-length. The calculations of grating intensity given in many 
text-books are so much waste paper. 

That the present theory is a ‘classical’ one, with dipoles radiating into 
a space obeying Maxwell’s equations, does not mean that the principles 
of quantum mechanics have been rejected. Quantum mechanics aims, 
in the domain of molar physics, at enunciating laws for the variation 
of statistical averages. It has not so far yielded the formulae we need 
here. Falling short of a theory of dielectrics, it has given a theory of 
dispersion in which an electron is acted on by its own atom, the incident 
radiation, and the reaction of the radiation. The constitutive equations 
of ponderable media in classical electrodynamics are built up by esti- 
mating the average effect of other atoms, and this will have to be done 
in the quantum theory before we can use it. 
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